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Abstract 

By constructing the Riemann surface controlling the resonance structure of Winter 
model, we determine the limitations of perturbation theory. We then derive explicit 
non-perturbative results for various observables in the weak-coupling regime, in which 
the model has an infinite tower of long-lived resonant states. The problem of con¬ 
structing proper initial wavefunctions coupled to single excitations of the model is 
also treated within perturbative and non-perturbative methods. 
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1 Introduction 


The only general analytic tool available up to now to study “realistic” quantum field theo¬ 
ries, such as for example Quantum Chro mo dynamics (QCD) in four dimensions, is pertur¬ 
bation theory. The great variety of the electromagnetic and hadronic phenomena observed 
at different energies may suggest that exact solutions will be beyond human capabilities 
for a long time, even though there has been some recent progress in QCD in the 1/Nc 
expansion, with Nc the number of colors [1]. In general, perturbative cross sections can be 
written schematically as: 

OO 

n=l 

where the cr„’s are real coefficients and g is the coupling of the model. QCD for instance, 
because of asymptotic freedom, is weakly coupled in the ultraviolet, where perturbation the¬ 
ory is therefore a natural tool, while it is strongly coupled in the infrared. The hadron mass 
spectrum, the scattering lengths, the string tension, the chiral symmetry breaking scale, 
the parton distribution functions, the power-corrections to high-energy cross sections and 
event-shape distributions, etc., are all well-known examples of significant physical quantities 
which fall outside the reach of perturbation theory. Furthermore, being unable to exactly 
evaluate the Un’s for arbitrary n, the problem of non-perturbative effects to observables is 
necessarily treated in a rather indirect way. In the last decades, many different approaches 
have been developed for this task: methods based on the summation of specific classes of 
diagrams (the 1/Nc expansion cited above, the renormalon calculus supplemented by the 
large-/^o limit, etc.), toy models in lower space-time dimensions (typically two or even three), 
numerical Monte-Carlo computations of the euclidean theory regularized on a lattice, direct 
comparison of perturbative cross sections with experimental data, and so on. We may say 
that a large part of high-energy theoretical activity of, let’s say, the last forty years has been 
devoted to gain some control on the non-perturbative effects. The problem is a recurrent 
one in QCD and it might have been a problem also in the (standard) electroweak theory in 
the case of a very heavy Higgs, let’s say mn ^ 800 GeV. In the latter case the scalar sector 
would become indeed strongly coupled. However, the Higgs boson was discovered at the 
Large Hadron Collider (LHC) in 2012 with a small mass, tuh — 125 GeV, as previously 
indicated by indirect measurements, so the problem of non-perturbative electroweak effects 
has at present a limited phenomenological relevance. In general, by using a (truncated) 
perturbative computation of a quantum field theory model to describe some specific phe¬ 
nomenology ranging, let’s say, from high-Tc superconductivity to strong interactions, one is 
always faced with the problem of ’’what is missing”, i.e. which effects he in the unevaluated 
terms or are actually ’’invisible” to perturbation methods. In such a situation, it may be 
interesting to study a model which, though not a quantum field, can be analyzed both in 
perturbative and non-perturbative way. It is in this spirit that we present a systematic 
study of the so-called Winter model [211H El E], a non-relativistic quantum mechanics 
model possessing, in the weak-coupling regime, an infinite tower of resonant states coupled 
to a continuum [7], with Hamiltonian in proper units: 


H 


dx‘^ 


H- 5{x 

Tig 


7l) 


( 1 ) 


2 



on the half line a; > 0 with vanishing bonndary conditions in the origin, 'ijj{x = 0,t) = 0. 

G M is the only conpling of the model. Winter’s model has an nnstable energy spectrnm 
in the free limit g ^ 0: for g ^ 0'^ the spectrnm is uniformly bounded from below by zero, 
while for g ^ 0~ the spectrum is unbounded because of the appearance of an eigenfunction 
in the discrete spectrum with energy e{g) ~ —l/(47r^5f^) —)■ —cxo (see figH]). The free limit 
is therefore a singular one. In classical terms, the particle falls in the potential trap located 
at X = TT. The instability is related to the behavior just of the fundamental state (i.e. of 
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Figure 1: Energy spectrum of the Winter model (arbitrary units). Repulsive case (g > 0) 
on the left of the vertical line (in red); attractive case (g <Q) on the right (in blue). Note 
the discrete line in the latter case. 


the state with the lowest energy) for g ^ 0, i.e. to the non-trivial “vacuum structure” 
in a small neighborhood of the free theory (see hgl2]). Furthermore, transition amplitudes 
for — 1 < < 0 contain contributions of the form coming from the bound state, non- 

analytic in the origin, which are typical of non-perturbative quantum field theory effects. 
A similar instability in Quantum Electrodynamics (QED) was conjectured in the 50’s by F. 
Dyson [8] 0. For 0 < a <C 1, where a = e^/(dvr) is the fine structure constant of QED0 with 
e the electron charge, the fundamental state is the vacuum, i.e. the ’’empty” state, without 
any electron-positron pair or any photon. Because of field fluctuations, electron-positron 
pairs (as well as photons) come out of the vacuum as virtual particles only, as their creation 
as real particles would increase the energy. On the other hand, for a < 0 particles with 
equal charges attract each other, while electrons and positrons repel each other, so that 
the Coulomb energy associated with an e’''e“ pair is positive. The rest energy of a system 

^ An instability mechanism similar to the QED one also occurs in the scalar theory, as well as in 
the simple case of a (quantum) anharmonic oscillator [^. 

^ Experimentally a = 1/137 at low energies. 


3 







Bound-State wavefunction 

fd(x, g) 



Figure 2: Wavefunction 'ipd{x,g) of the bound state of the Winter model for g = —0.07. 
It is also schematically shown the Dirac potential in x = n by means of the dashed vertical 
line. 


containing N pairs is 2Nmc^, with m the electron mass — proportional to iV — while 
the potential energy is of order |a|A^^/Ac — proportional to — where Ac = h/{mec) is 
the electron Compton wavelength. That implies that the potential energy overcomes the 
rest energy for large enough N. As a consequence, a huge quantity of e’''e“ pairs can be 
created as real particles out of the vacuum, which would then decay into a state full of 
electrons and positrons, by lowering its energy to arbitrarily negative values. The physical 
picture is the following: the e+e“ pairs are created close to each other, let’s say within their 
Compton wavelength (the QED interaction is local) and then, to minimize the energy, all 
the electrons fly on one side, coming close to each other, with the positrons flying on the 
opposite side. This spontaneous polarization of the vacuum should also occur in the usual 
particle states of the Fock spacqj. This instability of QED in a neighborhood of a = 0 
was related by Dyson to the divergence of the perturbative series in a. The spectrum 
instability of Winter model in a neighborhood of the free theory might suggest that also 
its perturbative expansion around = 0 is divergent — perhaps asymptotic to the exact 
theory for —)■ 0, as is supposed to be in QED. The reality is actually more intricate; after 
all, physical intuition is based on g eM., while a full understanding of the model, as we are 
going to show in detail, requires to complexify g. It turns out that the resonances of the 
model, which behave in space-time for x > 7i (i.e. outside the cavity) as 


g ikx—ik^t 

with k E C, are controlled by a multivalued function 

w = h{z) 

which is the inverse of the entire function 

e"' - 1 

^ = -, 

w 

^ The above argument is qualitative and does not provide any estimate of the decay time. 


( 2 ) 

(3) 

(4) 
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where w = 2Tiik and z = —g% The transcendental (inhnite order) mnltivalnedness of h 
is related to the fact that the model has an infinite tower of resonances. Each resonance, 
let’s say the n-th one with n a non-zero integer, is associated to a sheet S'„ of the Riemann 
snrface S of h; there is also an additional sheet, S'o, related to the bonnd state. Each S'^, 
with n 7 ^ 0 , has order-one (sqnare-root) branch points in 


Cn 


i 

27171 


(5) 


and at infinity, connecting Sn to S'o. That implies that S'o has a conntable set of order-one 
branch points {cnjnj^o accnmnlating at the origin as 


c„ —>■ 0 for n —>■ ±cxo. 


( 6 ) 


It also follows that different sheets, Sn and Sk with n, fc 7 ^ 0, ’’talk to each other” only 
indirectly, throngh S'o, which is a different (and more complicated) sheet with respect to 
all the other ones. Let ns also observe that the sqnare-root branch point is, in some sense, 
the most general conpling between different sheets one could think of. The perturbative 
expansions (series in powers of g) of quantities related to the n-th resonance (wave-vector, 
frequency, width, etc.) are convergent within a disk centered in = 0 of radius 


R 


n 


1 

27rn 


(7) 


The convergence radius is then controlled by the branch point c„ connecting S'„ to S'o. 
A non-zero radius of convergence for the expansion around zero was expected on phys¬ 
ical ground but, remaining in the physical domain G M, one could not have derived 
in a natural way its value. The instability of the spectrum, discussed above on physical 
ground, manifests itself mathematically in the fact that, according to eq.(| 6 ]), the origin of 
the “bound-state sheet” S'o is a non-isolated singularity. The above property is related to 
the specihc geometric structure of S': if the square-root branch points had coupled, for ex¬ 
ample, Sn to Sn+i for any integer n, the point g = 0 would have been instead an analyticity 
point for S'o. For quantities related to the bound state (such as its energy, its wavefunction, 
etc.) no convergent power-series expansion around g = 0 exists. While in the case of QED 
the instability of the vacuum is expected to manifest also in any state of the Fock space, 
in the Winter model it is restricted to the fundamental state. The second implication of 
eq.(l 6 ]) is that 

Rn ^ 0 for n —)■ ±cxd. (8) 

The consequence is that perturbation theory can accurately describe the dynamical prop¬ 
erties involving a finite number of resonances for small enough coupling, but it cannot 
describe quantities involving an infinite number of them. We will see that also the second 
class of observables contains fundamental physical quantities. 

Our motivation to further investigate Winter model might suggest that it is just a toy 
model for quantum held theory: that is not actually the case [lO]. Since it describes par¬ 
ticles conhned by potential barriers, with tunable efficiency, it is still currently used in 

^ The minus sign in front of g is inserted just for practical convenience. 
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quantum chemistry, together with its natural generalizations miDSlDS]. Let’s schemat¬ 
ically summarize the history of Winter model relevant to the present work. As far as we 
know, this model was originally introduced in [2], where the resonance properties of the 
spectrum for a small positive coupling were analyzed. The temporal evolution of metastable 
states of sinusoidal shape concentrated at t = 0 inside the cavity — the segment [0,7r] — 
was studied in [3]. In this work, post-exponential power corrections in time were explicitly 
calculated, conhrming a previous general analyticity argument about the breakdown of the 
exponential behavior at very large times. In ref. [5] Winter’s computation was repeated, 
hnding additional ’’non-diagonal” contributions to the exponential time evolution, which 
had been overlooked in |3]. These new terms imply a coupling of the initial state with all 
the resonances of the model and not just with a single one. The off-diagonal terms a small 
coupling 0{g) compared to the one in [5], but decay in general slower in time, dominat¬ 
ing then the wavefunction in a large temporal region. Actually, most of the non-trivial 
properties of Winter model — resonance mixing in particular — are consequences of this 
additional contributions to time evolution recently discovered. In [5] it was also found that 
the coupling of the initial state to the resonances was controlled at hrst order in g by an 
inhnite matrix of the form 

U{g) = l + gA, (9) 

with A a real antisymmetric matrix (see eq. fll32p L U{g) is then an inhnitesimal rotation in 
the inhnite-dimensional vector space of the resonances. It was then natural to conjecture 
that higher-orders in g would have led to the exponentiated form, i.e. to the unitary matrix 

U{g)=exp{gA) = l + gA+^A^+ ■■■ . (10) 

In order to check this structure, the second-order computation in g of U{g) was made in [6]. 
In addition to the expected term g‘^A‘^/2, it was also found a ’’large” term not compatible 
with any generalized form of exponentiation, whose physical interpretation was problematic. 
A possible investigation at that point could have been to push the perturbative expansion 
to 0{g^), in order to have some hint of the general structure and eventually resum the 
expansion at all orders in g, in the spirit of classical quantum-held-theory investigations. In 
this work however we follow a different route: we abandon perturbation theory and perform 
an analytic study in the complex plane of the functions controlling the expansion above 
for U{g). As we are going to show, that allows us to determine the convergence region 
of the perturbative expansion and to hnd explicit non-perturbative formulas for relevant 
observables which replace the perturbative ones outside their convergence region. 

This work is at the border of different areas: 1) quantum held theory, in particular high- 
energy theoretical physics, where the interplay between perturbative and non-perturbative 
effects is crucial, as already discussed; 2) mathematical physics, as we investigate the an¬ 
alytic and geometric structure of the Winter model and hnally 3) quantum physics — 
quantum chemistry in particular — where generalized Winter models are currently in¬ 
vestigated. Therefore, since the paper is of potential interest to readers with different 
backgrounds, we tried to be as simple and explicit as we could. The paper is organized as 
follows. In secj2] we summarize the properties of the spectrum of the Winter model and 
we discuss in particular the free limit —?• 0, in which the system decomposes in two non¬ 
interacting subsystems, a particle in the box [0, vr] and a particle in the half-line [vr, oo). All 


6 




that is in complete agreement with physical intuition. In secj^lwe schematically discuss the 
time evolution of wavefunctions initially concentrated inside the cavity and of sinusoidal 
shape. We also treat, in general terms, perhaps the most significant effect in the evolution 
of metastable states of Winter model: the above mentioned mixing of the resonances. In 
seclUwe specify the discussion on resonance mixing by using explicit power expansions in 
g and we discuss the problems associated with the perturbative expansion. In secEl we 
abandon perturbation theory and investigate the structure of the Riemann surface S of the 
multivalued function w = h{z) controlling the behavior of the resonances and the form of 
the inhnite mixing matrix U{g). In secj6] we re-analyze the properties of the resonances 
and of the mixing matrix by means of the exact (non-perturbative) information obtained 
with the previous geometric study. We deal again, in particular, with the problems encoun¬ 
tered with perturbation methods. Finally, in secJ7]we draw our conclusions and we discuss 
some natural developments of our work. There are also three appendices. In appendix 
we discuss the connection of our function h with the Lambert W function [T3], a kind of 
generalization of the complex logarithm. In appendix [B] we present a general treatment, as 
well as explicit formulas, for the branch points of the function h, which are crucial in our 
analysis. Finally, in appendix Owe show that the expansion for large \z\ of h{z), which 
involves the composition of the complex logarithm with itself, actually has exactly the same 
multivaluedness as h, which is just logarithmic, to 


2 Winter Model 


In general, the Hamiltonian operator of the Winter model reads: 

h? ( 9 ^ 

H = + XS{x-L), (11) 

2m ox^ 

where m is the mass of the particle, A is a (real) coupling constant and S{x — L) is the 
Dirac h-function with support in x = L > 0 [15]. The domain is the half-line 0 < x < cxo 
and we assume vanishing boundary conditions at zero0 

V’(x = 0,t)=0, t G M. (12) 


Formulas can be simplihed by going to a proper adimensional coordinate via 


and rescaling the Hamiltonian as: 



H = ^^H'. 
2mL^ 


(13) 


(14) 


The new (adimensional) Hamiltonian then takes the form in which it appears in the intro¬ 
duction: 

ox' ^ ng 

^ Equivalently, one may think to the problem in the whole real axis, x S R, with an additional infinite 
potential for x < 0. 
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and contains the single (real) parameter 


^ 2mXL 

(16) 

The time-dependent Schrodinger equation 


ih^ = Htl^ 
dt ^ 

(17) 

now reads 

dib 

(18) 

where time is rescaled as 

2mL‘^^ 

(19) 


Let us omit primes from now on for simplicity’s sake. It is possible to rescale the Winter 
Hamiltonian flTT|) in different ways, as made for example by Winter itself in [3]. The 
main point however is that we deal in any case with a one-parameter model describing the 
coupling of a cavity (the segment 0 < a; < tt) with the outside (the half-line x > vr). As we 
are going explicitly to show in the next sections and as anticipated in the introduction, for 
<C 1 resonant long-lived states inside the cavity come into play. 

2.1 Spectrum 

For g > 0 there is only a continuous spectrum, with eigenfunctions of the form [21 |3l |5] (see 
figsJS] and H]) 


V’(x; k,g) 


TT 


TT 4:a{k, g)b{k, g) 

-|- 9{x — tt) [a(/c, g) 


x) sin(fcx) -f- 

+ b(k,g)e-“‘’]\ 


( 20 ) 


and energies 


e{k) = /c^ > 0. 


( 21 ) 


The function 9{ii) = 1 for ?/ > 0 and zero otherwise is the Heaviside step function and the 
square root in eq. (l20|) is the arithmetical one, as a{k,g)b{k,g) > 0 for any k,g ^ W (an 
over-all phase is in any case irrelevant). The coefficients entering the eigenfunctions have 
the following explicit expressions: 


a{k,g) 

b{k,g) 



^ |exp(-2.«) - 1] ; 
[exp(+2.ii) - 1] . 


( 22 ) 

(23) 


These coefficients have the following two symmetries: 


a{-k,g) = -b{k,g); a{k,g) = b{k,g), 
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(24) 









Eigenfunction - Resonant Case 

k, g) 



Figure 3: Eigenfunction %l){x]k,g) for resonant k = I — g. Repulsive case (g = 0.05/- 
continuous (red) line; Attractive case (g = —O.OSj; dashed (blue) line. The enhancement of 
the amplitude inside the cavity is clearly visible, as well as the Oigpi^x = T^)/g) discontinuity 
of fj'^x = tt). 


where the bar denotes complex conjugation. The hrst equation says that '4){x-, k,g) is an 
odd function of k and implies that the zeroes of a are opposite to those ones of b, i.e. that 
if 

b{ko,g) = 0, (25) 


then 


a(-k„,g) = 0. 


(26) 


The second equation implies that the zeroes of a are the complex conjugates of the zeroes 
of b for conjugate coupling, i.e. that if 


b{ko,g) = 0, (27) 

then 

a {ko,g) = 0. (28) 

For real g, the zeroes of a are then the complex conjugates of the zeroes of b. From the 
two properties above it is possible to reconstruct for real g all the zeroes of a and of b, for 

example, from the zeroes of b in the forth quadrant. Finally, note that for real k and g 

a{k,g) = b{k,g), k, g eR, (29) 

so that 

a{k,g)b{k,g) = \a{k,g)\^ = \b{k,g)\‘^, k, g eR. (30) 

In eq. (l20|) we have assumed the standard continuum normalization: 

tpl^x; k', g) fj^x] k, g) dx = 6{k — k '), (31) 
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Eigenfunction - Non-Resonant Case 

lAC. !<■, «) 



Figure 4: Eigenfunction 'il!{x;k,g) for non-resonant k = 1.5. Repulsive case (g = 0.05/- 
continuous (red) line; Attractive case (g = —0.05/- dashed (blue) line. The amplitude inside 
the cavity is much smaller than outside it. 


with 5{q) the Dirac ^-function. The quantity k is a. real quantum number but, since the 
eigenfnnctions are odd fnnctions of fc, one can assnme /c > C@. Note that the spectrnm is 
bounded from below by zero, nniformly in > 0. As we are going to show in the next section 
and as anticipated in the introduction, that is no more the case for g < 0, because of the 
appearance of a discrete spectrum. The eigenfunctions can also be written in trigonometric 
form as [^: 


'ip{x;k,g) 




x) A{k, g) sin 


\kx\ + 9{x — tt) 


sm 


[kx + ip(k, 9 )]|, 


(32) 


The inside amplitnde is given by (see £g.(l5])) 


^ = , =. 
2-\/a(fc, g)b{k, g) a/ 1 + l/(7r5ffc) sin2/c7r + l/{27i‘^g‘^k‘^){l — cos2k7T) 

Note that0 

lim A{k,g) = 1, 

/c—>-oo 


(33) 

(34) 


as expected on physical gronnd: high-energy states do not see the barrier. The phase shift 
between the ontside amplitnde and the inside one reads (see hgE]): 


cot [^{k,g)] 


Tigk 

sin^(7rfc) 


cot(7r/c). 


(35) 


For the cotangent to be invertible, let us assume for its argument the following range: 

® The case fc = 0 has to be discarded as one obtains in this case, because of the boundary condition, the 
zero function, which is not an acceptable wavefunction. 

^ It may be argued that A{k,g), as a function of k, resembles the shape of the total cross section 
of electron-positron annihilation into hadrons as a function of the c.o.m. energy ^/s, a{e~^e~ —>■ h's), 
above a heavy quark-antiquark pair threshold. Roughly speaking, the quarkonium states correspond to the 
resonances of the particle inside the cavity. 
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Aik, g) 


Inside Amplitude 



Figure 5: Amplitude A{k,g) of the eigenfunction 'ip{x;k,g) inside the cavity. Repulsive 
case (g = O.lj; continuous (red) line; Attractive case (g = —O.lj; dashed (blue) line. 
Peaks become less marked with increasing order. The dotted horizontal line represents the 
asymptotic value of the amplitude at large energies. 


— TT < ip{k,g) <0 for > 0; 0 < ip{k,g) < vr for g < 0. (36) 

With the above choice of the argument, 

lim Lp{k,g) = 0, (37) 

fc—>- + 00 

i.e. there is no phase shift in the high-energy limit, in agreement with physical intuition. 

2.2 Strong-Coupling Regime 

Even though we are mostly interested to the weak-coupling regime -C 1 — to be more 
specihc, to subtle weak-coupling properties — let us briefly consider the strong coupling 
regime 

\9\ » 1- (38) 

Waves hnd in this case a small barrier at the point x = vr; in the strong coupling limit 
g —?• ±cx3, the Dirac potential completely disappears and the Winter Hamiltonian becomes 
then the Hamiltonian of a free particle, 

»o = -^. (39) 

on the positive axis, x > 0, with vanishing boundary conditions at the origin. Note also 
that 

A{k,g)^l and (p{k,g)^0 for g ^ ±oo. (40) 

A duality therefore exists between the strong-coupling regime of the Winter model S> 1 
and the weak-coupling regime of a particle on the positive axis subjected to a small potential. 
Let us remark that there is continuity in going from the strong coupling regime, ^ 1, 
to the weak-coupling one, -C 1. As we are going to show in the next section, the only 
singular limit is —)■ 0. 
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Vik, g) 


Phase Shift 



Figure 6 : Phase shift ip{k,g) of the eigenfunction 'f{x;k,g). Free case (g = 0): dotted 
(black) line; repulsive case (g = 0.035j; continuous (red) line; attractive case (g = —0.035/- 
dashed (blue) line. The curve for g > 0 has been shifted upwards by tt for comparison with 
the case g < 0 at small k. 


2.3 Weak-Coupling Regime and Free Limit 

We are interested to the weak-coupling regime 0 < 5 ^ -C 1 and to the free limit of the 
model S' —)■ 0^ [5]. It is clear that, unlike the previous case, taking this limit directly in the 
Hamiltonian, which has a pole in s = 0, is meaningles^. We therefore take the free limit 
on the eigenfunctions (which is instead meaningful) and then find which Hamiltonian has 
the limiting eigenfunctions. 

For k not close to an integer within 0{g), namely 

k^n + 0{g), (42) 


with n a non-zero integer. 


Hk,g)\ = \b{k,g)\ = O , (43) 

implying that the eigenfunctions have a small amplitude O {g) 1 inside the cavity — 

outside the cavity the amplitude is always 0{1), no matter which values are chosen for k 
and g, because of continuum normalization. With the generic values of k given by eq. fl42p . 
outside waves are not able to excite appreciably the cavity. In the limit s O’*" the 
eigenfunctions exactly vanish inside the cavity. 

® This is to be compared with quantum field theory, in which the coupling to be sent to zero usually 
appears in the numerator. The QED Lagrangian density, for example, reads 

C = - + b - m) V' + (41) 

where F^^, is the field tensor and the gauge potential of the electromagnetic field, is the Dirac field, 
e the electron charge, and the free limit is e —>■ 0. 


12 
















Let us then consider the phase behavior of the eigenfunctions. Unlike the discussion on 
the amplitude, let us hrst consider the free case. For g = 0, one obtains: 

(f = (ps{k, g = 0) = —Tik + STi, (44) 


where s is any integer. The eigenfunctions therefore read: 


k, g = 0) 


{-iy9{x 



k>0, kyn + 0{g), (45) 


with n a positive integer. In order to satisfy eq. (l36|) . the integer s has to be a function of 
k but, since an overall phase is not observable, one can take once and for all for example 
s = 0 , obtaining the usual eigenfunctions of a particle conhned to the half-line x > vr, with 
vanishing boundary conditions. 


'0(x = tt; fc, = 0) = 0. (46) 

The phase ip{k,g) does not provide any measurable information in the free case: it does 
not transfer any information from inside the cavity to outside it. In other words, the cavity 
is impermeable also as far as the phase (p is concerned. Let us now consider the interacting 
case, 0 < 1. For the non-exceptional values of k in eq. fl42D . it holds 

sin( 7 rA;) = 0{1) (47) 


and the additional term in ip{k,g) due to the interaction is just a small correction with 
respect to the free case: 


71 gk 

sin^(7rfc) 


0(9). 


(48) 


Therefore there is not any qualitative change of the prohle of (p(fc; g) in the region of generic 
/c’s given by eq. flT 2 l) (see hgj 6 ]). 

Let us now see what happens for 0 < (7 -C 1 and for (7 —)■ O'*" to the amplitudes and 
phases of the eigenfunctions for the complementary values of fc, i.e. for the “exceptional” 
values 

A: ~ n — ng, (49) 


with n any non-zero integer. Since 


\/l -I- TT^T?^ 

9),9]\ = \b[n{l - 9 ), 9] \ = -^-lf/| + O (/) , (50) 

the amplitude of '?/'(x; k,g) inside the cavity shows marked peaks. In the limit 5 ^ —)■ O’*" 

|a[n(l= |6[n(l-^),^]| ^ 0 (51) 

and the inside amplitude diverges. This divergence, which (at face value) is physically 
meaningless, actually signals a qualitative change of the spectrum. To have a hnite ampli¬ 
tude inside the cavity for g ^ 0, one has to impose that the eigenfunctions exactly vanish 
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outside it, necessarily obtaining states with a finite normalization. Therefore a discrete 
spectrum emerges out of the continuum one for the exceptional fc-values in eq. fHOD when 

—)■ 0 so that k ~ n(l — g) ^ n. (52) 


By normalizing these eigenfunctions to one, we obtain: 



x) sin(?7,x). 


(53) 


with n a positive integer. The above ones are the eigenfunctions of a particle confined to 
the segment 0 < a: < vr. 

As far as the phase behavior is concerned, the term on the l.h.s. of eq. fldSp has a double 
pole for k —>■ integer, while the independent term in eq. fl35p . cot( 7 r/c), only has a simple 
pole. The consequence is that the term proportional to g is not uniformly small in k even 
for | 5 f| -C 1 and actually dominates for k —)■ integer. In the region specified by eq. fl49p 

^ (54) 

for S' > 0 and g respectively. Furthermore, ip{k,g) passes through =F 7 r /2 with a high 
slope when k passes through n(l — g): 


*<==..(!-») i + \s) 


(55) 


This sudden phase variation is a typical resonant behavior [16]. An important point is that, 
unlike the free case, (p{k;g 7 ^ 0) is a continuous function of k and is a measurable. In the 
limit S' 0 the phase ip{k,g) develops infinite slopes at the integers k = n (see figj 6 |). 

By looking at eqs. fl4^ . fl50P and fl55l) . one may observe that the “effective coupling” of 
the n-th resonance to the continuum rather than g is actually ~ ng. That implies that 
one has to face strong-coupling phenomena for large n even when l^fl 1 , because the 
perturbative expansion involves powers of ng rather than powers of g\ we will make these 
considerations more precise in later sections. 

We can summarize the above findings by saying that, in the free limit, the system 
described by the Winter Hamiltonian decomposes into two non-interacting subsystems. 
The first subsystem is particle in a box, with Hamiltonian 


- 

Hi = 0 < Xi < TT, 

C/OC-^ 

with vanishing boundary conditions: 

■^1(2:1 = 0 , t) = '^1(0:1 = TT, t) = 0 , t G M. 
As well known, this system has an infinite tower of discrete states. 





sin (nxi ), 


(56) 


(57) 


(58) 
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with n G M+, and no continnnm spectrnm. The second snbsystem is a particle in a half-line, 
with Hamiltonian 


Ho = 




dxo ’ 


71 < X2 < OO, 


(59) 


with wavefnnctions vanishing at the only bonndary point x = it■. 


= TT, t) = 0, t G M. 


(60) 


H 2 has a continnous spectrnm only, with eigenfnnctions (normalized to a 5-fnnction) of the 
form 


(3)2; k ) 


y ^ sin [k{x 2 - vr)], 


(61) 


with k > Q. Therefore Winter Hamiltonian H{g) factorizes in the limit ^ 0+ in the snm 
of the Hamiltonian of a particle in a box and the Hamiltonian of particle in a half-line: 


lim H{g) =Hi + H 2 . (62) 

g^0+ 


2.4 Resonances 


To stndy the decay, as well as the formation, of metastable states, it is convenient to intro- 
dnce generalized eigenfnnctions Tp^x] k,g) with k E C, called resonances or antiresonances, 
which satisfy pnrely ontgoing or pnrely incoming bonndary conditions, i.e. 

'ip{x;k,g) = 0 for x > ti, (63) 

with 

Re A: > 0 or Re A: < 0 (64) 

respectively. Let ns remark that while in the case of the Winter model the bonndary 
condition above, implying 

7/j{x-, k,g) Tii (65) 

can be imposed at any point ontside the cavity, for a general short-range potential, one 
has to impose it at x —)■ -|-oo. By eqnating to zero the coefficient b{k,g) of the component 
exp{—ikx) in '0(a;; k,g) (see eq. fl20|) ). we obtain the transcendental eqnation in k 

exp(27rAA:) -|- g2nik —1 = 0, (66) 

having a conntable set of solntions {k^^\g)} lying in the lower half of the A:-plane, 

lmk^"‘\g) < 0 (67) 



and reading for \g\ -C 1: 

k^^^g) = n{l-g + g^) - iirn^g'^ + O {g^) , (68) 

where n is any non-zero integer. For n > 0, the zeroes lie in the forth qnadrant and are 
associated to resonances, while for n < 0 they he in the third qnadrant and are associated 
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to antiresonances. In general, the fnnctions with n 7 ^ 0 are defined as the zeroes 

of the transcendental eqnation above, 

b[k^-\g),g]^0, ( 69 ) 

satisfying the initial condition 

k(^\g = 0)=n. (70) 

These fnnctions are fnndamental elements of the Winter model. In the next section we will 
present a higher-order pertnrbative expansion for k^'^^ (g) — as already noted, we will find 
that the expansion parameter, rather than g, is actnally gn. In sec. ([5]), in order to fnlly 
nnderstand their properties, we will analytically continne the k^'^\gys for complex g. All 
the k^'^\g)^s will tnrn ont to be the varions branches of the mnltivalned fnnction k = k{g) 
implicitly defined by b [k{g),g] = 0 . 

Becanse of the first symmetry between the fnnctions a(fc; g) and b{k] g) discnssed in the 
previons section, the zeroes of fho fnnction a{k;g) have a positive imaginary part 

for any n 7 ^ 0. They can be defined as 

= -k^'^Kg) = -n{l - g + g^)+ iirn^g^ + ■■■ (71) 

and are associated to resonances for n > 0 (second qnadrant) and to antiresonances for 
n < 0 (first qnadrant). Ronghly speaking, becanse of the first symmetry between the 
coefficients, all the dynamical information is already contained in jnst one coefficient; we 
have therefore considered explicitly only the coefficient 6 (fc; g). By using also the second 
symmetry, one can easily show that for (7 G M the real part of k^'^\g) is odd in n, while the 
imaginary part is even in n. From the relation 

X"('j) = (72) 


it follows indeed 

{ 9 ) + { 9 ) = - {g) + {g) , ( 73 ) 

where k^'^\g) = k^\g) + ik^\g) with k^\g), k^\g) G M. For real g the Winter Hamilto¬ 
nian is real and therefore time-reversal invariant, so that the formation and the decay of a 
resonance occur in the same way. 

According to the formula above, for n = 0 one would obtain k^^^g) = 0; unlike the 
other cases, the function k^^^g) does not posses a convergent expansion in powers of g. 
Relevant expansions in this case are an expansion in powers of 5 ^ -|- 1 for + 1| -C 1 

? 0,0 

- -(9 + 1) + 5-(g + 1)^ + 0 [(g + 1)=] . (74) 

TT OTT 

as well as an expansion involving exponentials and poles for g < 0 , \g\ < 1 : 


k^^\g) 



T/9 


+ C> (e^/^) 


( 75 ) 


We will see the physical significance of these expansions in the next section. Let us remark 
that k^^^g) is a purely imaginary number for g G M“, with a positive imaginary part for 


- 1 < ^ < 0 . 


(76) 
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The latter case corresponds to an exponentially decaying wavefnnction for x —)■ +cxo with 
negative energy, i.e. a bonnd state. 

Coming back to the resonances, we obtain for their wavefnnctions (n > 0): 


and for the (complex) energies 

The resonances evolve diagonally in time by means of the factors 


g) ^ exp[-ze^-\g)t] 


= exp 




X 1 

X 

J 

(g) t] 

( 77 ) 


( 78 ) 


( 79 ) 


Since the energies are complex for 7 ^ 0 , on the last member we have split them into real 
and imaginary parts as: 

9 '”)( 9 ) = a.(">(9) - |r>'‘>(9), (80) 

where is the freqnency and is the decay width of the resonance n: 

a;W(( 7 ) = Re = n\l - 2g) + O (/) ; (81) 

r(’^)(c/) = -2lm[k^^\g)^] = 47igV + O {g^) . (82) 

Note that u^"'\g) is even in n, while r^"^( 5 f) is odd in n, and that = 0; ( 7 ) = 1. Let 

ns make a few remarks. 


1. The exponential decay of a resonance with time is not in contradiction with the 
conservation of probability (i.e. of the number of particles) because {x,t; g), as 
a function of x G M’*', is not a normalizable wavefnnction and describes an outgoing 
flux of particles at x — +cxd [^. Actually, a resonance wavefnnction is not even a 
bounded function — like the ordinary eigenfunctions— and diverges exponentially 
for X —)■ +CX 30 since Imk^'^^g) < 0 |^. By writing indeed 

k = ki — ik2, ki, k2 G M, (83) 

the outgoing boundary condition and the positivity of the width, 

fci > 0, r = 4:kik2 > 0, (84) 

imply /c 2 > 0 , so that the resonance wavefnnction diverges exponentially for x —)■ +00 
as 

ikx _ ikix-\-k2X. 

o o ^ 

® For real < 7 , the Winter Hamiltonian is hermitian only with boundary conditions implying no net flux 
of particles at infinity. 

In the large-time numerical evolution of wavepackets initially concentrated between potential wells, 
such an exponential increase is actually observed in a large space interval mm- 
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2. The resonances, just like the ordinary eigenfunctions, are smooth functions in M+\{ 7 r}, 
while they are only continuous at a: = vr. Even a finite discontinuity would indeed 
produce an infinite average kinetic energy (see next section); 

3. Outside the cavity, is 0{g) compared to the inside, because of the explicit 

g factor in the second term on the r.h.s. of eq.dTTIb As expected on physical ground, 
apart from the (slow for -C 1) exponential divergence for a; —)• +cxo discussed above, 
the wavefunction is concentrated inside the cavity. 


2.5 Attractive Case 

Even though we are mostly interested to the weakly repulsive case, i.e. to 0 < -C 1, let 
us briefly discuss the modifications of the spectrum in the attractive case, i.e. for (7 < 0 [5]. 
In the latter case, there is a continuous spectrum as in the repulsive case, while for 


- 1 < ^ < 0 


( 86 ) 


there is also a discrete spectrum consisting of a single bound state. The discrete (d) eigen¬ 
function reads: 


g) = C, [e(T! - X) - e-i's**) + e(x - it) - l) 


(87) 


and has the negative energy 


£d{.g) = -x^{.g) < 0 . ( 88 ) 

The quantity xig) ^ is the imaginary part of the bound-state solution of the equation 
b{k,g) = 0 , which is purely imaginary, 

k^°\g) = ix{g)- (89) 

It is the positive solution of the transcendental equation 

g-27rx(3) = 1 + 2Tigx{g) . (90) 


By normalizing the wavefunction to one, the constant above reads: 


C, = 


x{g) 


g 27 rx(g) — 1 — 2'n'x{g) 


(91) 


For a “loosely-bounded” particle, i.e. for |1 + S'] -C 1, the transcendental equation above 
has the approximate solution (see eq. fTMll ) 

x(g) - ha + 1) + |-(J +1)^ + 0 [(g + 1)3] , (92) 

71 071 

while for a “tightly-bounded” particle, i.e. for a negative coupling of small size, g < 0, 
\g\ 1 , one has the expansion (see eq. flTSjl l 


x{g) 


1 

27ig 


T/9 


+ C> 


(93) 
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We will rederive the above expansions when studying the analytic continuation of the 
functions k^'^\g) controlling the resonance behavior of Winter model (see later). In the 
case <C 1 (tight binding), by omitting exponentially small terms, we have the explicit 
formula: 


g) ~ 


1 N f X - 7r\ . f TT - X 

— , &(7r — x) exp ——- + U[x — tt) exp ——- 

\M9\) ' ’ V2>risi 


(94) 


and 

- -J^r (“) 

As discussed in the introduction, for g ^ the spectrum becomes unbounded from below, 
while it is uniformly bounded by zero for —)■ O’*". 

Let us now consider the free limit in the (more complicated) attractive case, 5 ^ —)■ 0“. 
Roughly speaking, in the this case, factorization is not so clean because of the presence of 
a discrete spectrum. As we have seen, for —)■ 0 the continuous-spectrum eigenfunctions 
leave the cavity for generic k values, while they concentrate inside it for the exceptional k 
values k ~ integer. On the contrary, the bound-state wavefunction extends symmetrically 
on both sides of the potential wall for g ^ 0~ (see hgj2]). However, since for g ^ 0~ the 
wavefunction is concentrated in a thin layer around a: = tt, of width ^ 47r|5f|, the resulting 
coupling between the box [0, tt] and the half-line [tt, 00 ) is small. In a weak sense indeed: 


while 


lim ijlix-^g) = (5(x - vr), 

g-)-0- 


lim tlJd{,x]g) = 0 . 


(96) 

(97) 


3 Temporal Evolution of Metastable States 


We study the forward time evolution, t > 0, of wavefunctions t\ g) which coincide at 

the initial time, t = 0 , with the box eigenfunctions and vanish outside it [31 El [ 6 ]: 


= 0 ; g) 


\/2ln sin (/ x) for 0 < a: < vr; 

0 for TT < a: < cx), 


(98) 


where I is a positive integer. The above wavefunction is just a wavepacket with support^ 
entirely contained inside the cavity. Physically, that corresponds to consider an initial 
state containing an unstable particle such as, for example, a Z^, but not any of its decay 
products. By means of a spectral representation in (ordinary) eigenfunctions and contour 
deformation in the /c-plane, the wave-function at time t > 0 can be exactly written as [5]@ 


i^^’'\x,t;g) = '^V{g)in0^^\x,t;g) + P^^\x,t; g), (99) 

n=l 

The support of a numerical function / : D —>■ C is the closure of the set where the function is not 
vanishing: Supp/ = {x € D ■. f{x) A 0}. 

For — 1 < g < 0 there is also an additional term on the r.h.s. of ea. IMl) coming from the discrete 
spectrum, a case which we do not consider explicitly. 
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where 6^^'>{x,t] g) is the n-th resonance wavefunction constructed in the previous section 
and g) is a non-exponential (power-like) contribution, having the exact integral 

representation 


P^^\x,t] g) = e < 6 *( 7 r — a:) j (^ke ^f) sin (fc e 


^dk + 


9{x — tt) 


qii) g) a{ke-^^/^; g) P)kx-kH 


.0 

00 


+ / gl'* (k ij) b(k g) dk 


( 100 ) 


with 


q^^\k] g) 


(— 1 )^/sin fcvr 1 

F -/2 4a(fc; g) b{k- g) 

(-l)^/sin kn _ 1 _ 

kP' — P 1 + 1/ {T^gk) sin 2kTT + 1/ {2TT'^g‘^k‘^){l — cos 2 /c 7 r) 


The dependence on the initial state, i.e. on /, is contained in the hrst factor on the second 
and last member of the above equation, while the dependence on the dynamics is contained 
in the second factor. The following asymptotic expansion for f S> 1 holds: 


P^^\x,t-, g) ^ 


eW4 ('_iy g 

\/2 I 1 + 17 



- 1 “ 9[x 





( 102 ) 

The following remarks are in order: 1) the above function is continuous in a; = tt, with a 
discontinuous first derivative at the same point; 2 ) the power corrections are O {g'^) inside 
the cavity and O {g) outside it, implying larger contributions in the latter case for -C 1 . 
Power contributions however vanish in both cases for ^ 0, while for g —)■ ±cx) they 
represent typical dispersive behavior. As we are going to show in the non-perturbative 
section, in the strong-coupling limit the resonance contributions indeed disappear, because 
the functions k^'^\g) have imaginary parts ^ — cx) [5]. 

In the following we shall not concentrate on the post-exponential power-corrections in 
time related to PP\x,t; g), which do not have a resonance interpretation and therefore are 
not relevant for the present discussion. Physically, they involve the emission of very low- 
energy particles out of the cavity at very large times. We will also take t sufficiently large 
to avoid the considerations of the pre-exponential effects. The latter are related to a fast 
rearrangement of the initial wavefunction, which modihes its short-wavelength components 
in order to “£t” inside the cavity [3]. For 0 < 5 ^ 1 there is a large temporal window 

where the exponential decay is a good approximation inside the cavity [SIE]: 


1 < t < 


log(l/^) 

9^ 


(103) 
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Finally, the quantity V{g) is an infinite matrix describing the coupling of the initial 
state to the resonances, with entries 


( . ^ (-1)^+^ 21 g exp [ink^^^g)] 

^ /c(»i )(^)2 _/2 exp [27iik^^'>{g)] + g 


(104) 


where in the first factor on the r.h.s. we have isolated the dependence on the initial state, 
i.e. on /. As we are going to explicitly show in the next sections, the matrix elements 
V{g)in, even for very small g, decay quite slowly for l,n ^ oo. However, at fixed x and 
f > 0, that does not produce any convergence problem on the r.h.s. of eq. fl^ . since the 
resonances g)^s decay exponentially with n. Their widths, T^^\g)^s, grow indeed 

faster with n than the imaginary part of the k^'^^gYs (see eqs. (12321) and fl235p h 

r<">(g) nlogn while lmk^^\g) ^ logn (n S> 1). (105) 

In other words, the series on the r.h.s. of eq. (l99D is convergent even for a slow decay of 
V{g)in for n —)■ oo as long as f > 0. 


3.1 Resonance Mixing 

Let us now discuss in general terms the content of eq. (11041) . which is the second funda¬ 
mental element of the Winter model. In the next sections we will be more concrete by 
inserting explicit approximations, holding in different regimes, for the k^"'\gys inside of 
V{g). Eq. (11041) is indeed an exact relation for the matrix entries of V{g) in terms of the 
coupling g and of the functions k^^\gys. Unlike the dependence on n, which involves the 
detailed form of k^'^Yg), the dependence on I, i.e. on the initial state, of V{g)i^n is quite 
explicit; in particular: 

V{g)i,n^j for / > n. (106) 

Since the functions k^^\g) are continuous in = 0 (being analytic), the diagonal elements 
V{g)n,n are enhanced for |5f| <C 1, because in this case k‘^^\g) ^ n and a small denominator 
arises. That is in agreement with physical intuition: one expects the l-th initial wavefunc- 
tion in eq. (l98|) to excite mostly the l-th resonance. A fundamental point is however that the 
above property is not exact: the initial box eigenfunctions do not excite a single resonance 
at a time, but in principle the whole spectrum, with strengths given by the elements of the 
matrix V{g). Since the k^'^YgYs, being analytic functions, cannot have non-isolated zeroes, 
it follows that the V{g)i^nS are all non-zero for most of the values of g. Realistic physical 
models have many excitation^ so, in order to have real control on the system, one is faced 
with the problem of constructing initial states which excite, as much as possible, a single 
state at a time. In experimental high-energy physics, for example, that is the problem of 
producing monochromatic beams of specified particles, such as antiprotons or neutrinos, out 
of collisions. In lattice QCD, that involves the construction of a good interpolating field, 
exciting out of the vacuum, within human capabilities, only the states one is interested 
to. In essence, the problem is that of eliminating, as far as we can, the backgrounds. Let 

Just think to the QCD description of the collisions at the LHC, where dozens of hadrons are created 
at each event! 
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US stress that without some background subtraction there is no point even in building up 
experiments. In our case, to excite one resonance at a time, an idea could be to construct a 
new initial wavefunction, i.e. a new wavepacket, by inverting the inhnite matrix V{g) and 
applying it to the box eigenfunctions. We then go from 


= 0 ; g) 



x) sm{lx) 


(107) 


to 

f2 °° 

<!>^^'>{x,t = 0] g) = J-e{7r-x)'^{V{g)-^)^^sm{nx). (108) 

^ ^ n=l 

That way one obtains a formal initial wavefunction evolving, as far as exponential evolution 
is concerned, exactly as a single resonance: 

^^^\x,t-,g) = e^^\x,t-,g)+ U^^'>{x,t-,g), f > 0, (109) 

where 

OO 

U<^^\x,t;g) ^ Y. 9)- (HO) 

n=l 

The “experimental meaning” of eq. fllOSp is rather transparent: in order to excite the l-th 
resonance only and then observe a ’’diagonal” time evolution as in the free case ( 5 ^ = 0 ), 
one should prepare the initial state as the coherent superposition of free eigenfunctions 
given at its r.h.s.. To dehne invertibility for an inhnite matrix (and eventually give a 
rule for explicitly computing the inverse) requires in general a delicate limiting procedure. 
We will deal with this issue in the next section by means of perturbation theory and we 
will give some results concerning the exact inverse in the last section of the paper. The 
main problems in the above strategy to identify the excitations of the Winter model, i.e. in 
eq. fllOSp . are however the following. As already discussed, 9^’‘'>{x, t] g) diverges exponentially 
for X —>■ +CX 0 , so it is not even a bounded function. 


^ (Ill) 

while the “counter-rotated” initial wavefunction f = 0 ; 5 ^), being a wavepacket with 

support inside the cavity, is supposed to be square summable, 

= 0;c/) G (M+) . (112) 

However exact time evolution preserves the norm of the states, so ^^^\x^t]g) should be 
square summable at all timeso 


^^’■\x,t-g) e (M+) . 


(113) 


On physical ground, it seems difficult to completely ”£11” a resonant state by evolving a 
properly constructed wavepacket. The physical idea behind the use of resonant states to 

In principle, non-exponential corrections could invalidate this reasoning but, at the present level of 
understanding, that does not seem to be the case. 
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describe the long-time evolution of specific wave-packets (i.e. initially concentrated between 
potential wells), is indeed to consider a limited region of space, even though the problem is 
formally defined on the whole positive half-line. On the technical side, the problem of the 
above strategy is that the matrix elements V{g)in, even for very small g, decay quite slowly 
for l,n ^ oo. By inverting the matrix V{g), one obtains an infinite matrix V{g)~^ with 
elements (y{g)~^)i^ also slowly decaying for /,n —)■ oo — at least, that is what happens in 
perturbation theory and what exact (numerical) computations point out (see next sections). 
When one applies the matrix V{g)~^ to the inhnite vector containing the initial data. 



sin X, sin 2 x, sinSx, 


). 


(114) 


singular initial wavefunctions are obtained. A final remark about the formal construction 
above is in order. Eq. flinsp represents the initial wavefunction f = 0;^') as a sine 

Fourier transform with coefficients {V{g)~^)in: the set {sinnx, n G M} forms indeed a 
complete orthogonal basis on the interval [ 0 , 7 r] [1^. It is clear that one could have taken 
any other basis, as the results are independent on that (arbitrary) choice. 

In view of the above facts and considerations, it is natural to take a more ’’modest” 
approach to identify the excitations of Winter model, consisting in ” counter-rotating” the 
low-energy resonances only. Let us then assume that one is interested to the first N reso¬ 
nances only and rewrite eq. fl99p in the form 


N 

(x, t;g) = '^W{g)in 9^^^ (x, t; g) + r(x, fi g,N) + (x, t; g ), (115) 

n=\ 

where 1 < / < and W{g) is an N x N matrix with entries equal to the corresponding 
ones in V (g): 

W{g)in = V{g)in, /,n = 1 , 2 , • • • , A^, (116) 

while 

OO 

g,N) = V{g)in9^^\xfi-g). (117) 

n=N+l 

The latter is a remainder function containing the high-energy resonances which one does not 
want to treat explicitly. Since the lifetimes of the resonances quickly decay with increasing 
energy, we expect that neglecting the remainder function on the r.h.s. of eq. flllSp should 
be a reasonable approximation for t sufficiently large — but of course still within the 
exponential time region specified by eq. (11031) . Furthermore, since n 7 ^ / on the r.h.s. of 
eq. flllTj) . r9fix,t; g,N) involves small 0{g) off-diagonal couplings of the resonances to the 
l-th state (see next section). The meaning of the truncation at N should be rather clear. 
For A^ = 1, for example, one explicitly considers the first resonance only; an example of this 
approximate scheme with N = 2 has been sketched in [ 6 ]. We then consider the evolution 
of the wavefunction 0*^^^(x, f; g,N) given at the initial time, t = 0, by: 


(j)9)(x,t = 0; g,N) 



9{71 — x) 


^)z„sin(nx). 


n=l 


(118) 
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Its evolution reads: 


g,N) 


e^^\x,t-,g) + g,N) + 7r®(a;,t; g,N), 


(119) 


where 


N 

p^^\x,t] g,N) = '^{W{g)-^)^y^\x,t-, g,N) 

n=l 

oo 

= Y. RU9,N)9^^\x,t-g), ( 120 ) 

s=N+l 


with 


and 


N 

n=l 


N 

n<'\x,t-,g,N) = j). 

n=l 


( 121 ) 


( 122 ) 


The function p^’‘\x,t; g,N) represents a contamination coming from the resonances with 
order greater than N, to the l-th resonance state {I < N). As already explained, at 
sufficiently large times, such a contamination should be negligible, because the high-energy 
resonances have small lifetimes and small couplings to (j)^’‘'^{x,t = 0; g,N). 

We will review the perturbative computation of the matrix V{g), as well as of its inverse, 
in the next section, and we will present some elements of the exact theory in the hnal section 
of the paper. 


4 Perturbative Analysis 


By perturbative computation we mean an expansion in powers of g, for | 5 f| -C 1, of all the 
relevant quantities: wavevectors, frequencies, decay widths, mixing matrices, etc. In order 
to simplify the formulas and identify some structure, if any, it is convenient to allow for an 
ad-hoc normalization of the resonances. That is accomplished by dehning the ’’renormalized 
resonances” out of the starting (or ’’bare”) ones: 




S‘">(x,i;g) 

Z(»)(j) 


(123) 


where Z^'^\g) is a function of n and g whose value will be specihed later. Eq. flQ^ is then 
rewritten as [6] 

CXD 

ij^^\x,t]g) = YRin{9)i^'"\x,t-g ), (124) 

n=\ 

where the renormalized mixing matrix U{g) is given, by consistency, by (we have just 
multiplied and divided ^^'^\x,t;g) by the symbol Z^^\g)): 


U{g\n = V{g\nZ^-\g) (no J^). (125) 

n 
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We then insert the small-^f expansion for the inside the functions entering the r.h.s. 

of eq. fl99|) : 

= n — ng + (^n — g^ + — n^g^ + O (^g^) . (126) 


For the frequencies and widths, we obtain: 

J^\g) = n^l-2g + 3g^) + O (g^) ; 

= 47m^/(l - Ag) + O . 

By choosing for instance [6] 

Z^^\g) = 1 + I + QiTrn -^g"^ + 0 {g^) , 


U(g) reads in matrix notation, up to second order included: 

U(gr) = 1 + 5'r-^ + + O [g ^), 

where we have introduced the renormalized coupling 


9r = g-^g'^ + 0 {g^) . 


A is a real antisymmetric matrix with entries 

'4'." " f” 

and H is the real diagonal matrix with elements 


/ ^ n. 


Ai,i = 0 


Hl,n 


^ ^l,ni 


(127) 

(128) 

(129) 

(130) 

(131) 

(132) 

(133) 


where 6i^n = 1 for I = n and zero otherwise is the Kronecker delta. Let us remark that since 
for large n 


sin \k'^'^\g)x\ 


_g7r3^7i,2a; 

2 


while ss 


(134) 


the series on the r.h.s. of eq. fIMll is convergent for any f > 0. In secj6] we will see that 
a similar convergence mechanism also holds in the non-perturbative case (schematically: 

—)■ Inn while n^ —)■ nlnn). 

As already discussed, it is not straightforward to dehne the inverse of an inhnite matrix, 
as it not a matter of linear algebra (the determinant, for example, is not defined). However, 
a formal inverse of U {gr) exists in perturbation theory and is given, again up to second 
order in gr, by: 

U~\gr) = 1- gr A + ^grA‘^ - inglAH + O {gfj . (135) 
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Let us then consider the initial state 


OO 


ip^^\x,t = 0; g) 


9{ti - x)J2iU{9) sin(na:). 

n=l 


(136) 


The initial wavefunction is naturally decomposed in three structurally different contribu¬ 
tions: 

ip^^\x,t = 0; g) = a^^\x; g) + g) + g). (137) 

The first contribution is obtained by omitting the last O {g"^) term on the r.h.s. of eg. (11351) 
and reads: 


a^^\x- g) ~ Y ^ (l - I) - x) sin [/(I - g) x] 

~ e{'K - x)i^^\x,t = ^] g), (138) 


i.e. it is basically the renormalized resonance at t = 0 inside the cavity, continued to zero 
outside. The physical interpretation of the above equation is quite appealing: if you want 
’’diagonal” time evolution, you have to prepare the initial state with the wavevector k^^\g) 
of the l-tli resonance, not of the box eigenfunction k^’'\0) = 1. The problem is that the 
above function has a finite 0{g) discontinuity at the border of the cavity, 

lim a^''\x; g) ~ (139) 

X—>-7r“ 

while the ordinary eigenfunctions or the generalized eigenfunctions (resonances) are con¬ 
tinuous, so there is a substantial loss of regularity. Actually, a discontinuity produces an 
inhnite average kinetic energy. Schematically, with 


we have 


where 0 
we have 


'ip{x) ~ 9{tt — x) for a; ~ TT, 


(140) 


indeed 


dtl) 


(T) J \^) ~ ~ 

< £ <C 1. With the natural regularization coming from the Fourier representation. 


m our case: 



dx ~ 2'Kl‘^g‘^ N, 


(142) 


where we have truncated the sine series to the A^-th harmonic; one could choose for example 
N ^ l/^f, obtaining then an 0{g) correction to the leading-order kinetic energy. Further¬ 
more, the average potential energy of the Winter Hamiltonian, which involves the operator 
5{x — 7i), is not well defined for a wavefunction with a hnite discontinuity at a; = 7rD£|. That 

Formally: 

g) e (M+) but g) ^ (K+) , (143) 

where (R'*’) is the Sobolev space of the square-integrable functions on R+ with square-integrable (weak) 
first derivative again in K+ [TH] . 
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however is not the whole story and an even greater concern is with the term —ing'^AH. 
The latter comes from the imaginary term in the expansion above for (g) and 

reads: 

(— 2ln? 

{AH)i ^ = - -—— - for / 7 ^ n and 0 otherwise . (144) 

Since the latter has a large size for n S> /, it is natnrally decomposed as: 

= (-1)"+'+! -Jh- + (-l)»+'+i 21 . (146) 

The second contribntion to the initial wavefnnction, coming from the hrst term on the r.h.s. 
reads, omitting the trivial factor ^2Itt 9{tt — x): 

g) = —2'Kigl? In (^cos ^ j sin(/a;) + • • • , (146) 

where by the dots we mean less singnlar terms for x ^ 7r~ with respect to the one explicitly 
written. Unlike the ’’renormalized contribntion” a^^\x;g), the above fnnction is continnons 
for X —?• 7 r“. Actnally, it is proportional to the lowest-order term sin(/x), with a proportion¬ 
ality constant diverging logarithmically for x ^ Tr~. The best interpretation of this term 
we could think of, is that of some kind of imaginary contribution to the renormalization 
constant Z^’'\g), diverging for x —)■ 7r“, i.e. approaching the border of the cavity from the 

insidJ^. 

The third contribution to the initial wavefnnction is the one coming from the second 
term on the r.h.s. of eq. (ll45p and is the most singular one: 

9) = tan (|) + ''' • (147) 

The function above is not proportional to the lowest-order wavefnnction sin(/x) and 
contains a power divergence ~ l/( 7 r — x) for x —)■ vr”. Let us remark that one cannot even 
interpret eg. 01471) in a weak sense. In the latter case one should indeed resum the involved 
series as: 

00 ^ 

^(_l)«^+i sin(? 7 ,x) = - P.U. tan —, (148) 

2 2 

where P.U. denotes the principal value. Since the expression on the r.h.s. has then to be 
multiplied by 9{7i — x), one ends up with an ill-dehned distribution product. Then 7 (d(x; g) 
is not locally integrable and, as a consequence, not square integrable. The conclusion is 
that 7 *^^^(x; g) is totally unacceptable from physics viewpoint. 

To summarize, without further input, the best one could do is to use the truncation 
method described in the previous section to identify operationally the excitations of the 
system. One could arbitrarily choose some N and then see the dependence of the results 
on the truncation order to hnd some stability interval for N, if any. Actually, in the next 
sections, by means of a global study of the function k = k{g), satisfying b[k{g),g] = 0 for 
any g E C, we will be able to determine the natural truncation order of the perturbative 

In general, an imaginary contribution to the (on-shell) renormalization constant Z of a quantum field 
signals instability of the related (massive) particle [TO] . 
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results and to compute some non-perturbative effects as well. The conclusion we will reach 
is that for fixed g 0 the above formulas turn out to be applicable for small N only: 


N < 


1 

27r|^r 


(149) 


For larger N, one has to use quite different formulas for the mixing matrix V{g) (or U{g)), 
which have a completely different asymptotic behavior. 


5 The Multivalued Function h{z) 

By defining 

z = —g] w = 2m k, (150) 

the problem is, as discussed in the previous section, that of inverting the entire function 

/ : Cw —t Cz (151) 

dehned, according to eq. flOOl) . by 

^ = /(w) = ^-- (152) 

w 

for ta 7 ^ 0 with /(O) = 1. We want indeed to know the possible wavevectors A: as a function 
of the given coupling g, i.e. at hxed physics. Since f{27iin) = 0 for n ^ 0 integer, / is not 
injective (one-to-one) and the formal inverse h is not single-valued. We then proceed to its 
analysis by means of the following steps: 

1. We identify the branch points of h by means of a local study of the inverse / [20] : 

2. We perform local expansions of h at various points and we match them by evaluation 
on the overlaps of the convergence regions (whenever not empty), or by means of 
proper paths connecting them; 

3. We cut the complex ^-plane and then we glue the sheets in order to hnally construct 
the Riemann surface of h. 


5.1 Branch Points 

When the hrst derivative of z = f{w) vanishes, let’s say in w = d*, this function is not 
one-to-one in a neighborhood of d*, so the inverse 

h:Cz^ C^, (153) 


w = h{z), is multivalued in a neighborhood of c* = /(d*). The derivative 





vanishes when w 7 ^ 0 is a solntion of the transcendental eqnation 


= l-w. 


( 155 ) 


By nsing eq. fll52p . the above eqnation can be transformed into the following eqnation in 

(156) 




The order of the branch point can be determined by evalnating the higher-order derivatives 
of f{w) in d*. It holds 


f'M = 


{w^ -2w + 2)e^ -2 




so that 


rHifH=o- ^0. 

Since the second derivative does not vanish, we have that for |t(; — d*| -C 1 
z — ^ {w - ^ w - K, {z - 


(157) 

(158) 

(159) 


SO w = h{z) has a hrst-order branch point in z^. Let ns call {cn}nj^o the set of the solntions 
to eq. fll56l) different from z = 1 (which is not a branch point). The relation between the 
images of the branch points and the branch points c„ = / (d„) reads: 

Cn = ^ ^ , , n 7 ^ 0 . (160) 

1 (2j2 

As shown in appendix [HI all the c„’s lie in the hrst and the fonrth qnadrant of the z- 
plane (see hg|7]). Let {cn}n>o be the solntions on the hrst qnadrant, ordered according to 
decreasing distance from the origin. 


|ci| > |C 2 | > [cal >.> |c„_i| > |c„| > \cn+i\ > • • • > 0, Vn > 1. (161) 

Since eq. (11561) has real coefficients, solntions come in pairs of complex-conjngated points 


Table 1: Exact and asymptotic values of the first few branch points Cn’s (of order one) of 
h, together with the convergence radii Rn’s of the power expansions of the hn{z) ’s around 
z = 0 and the convergence radii pn’s of the semi-integer power expansions around the Cn’s 
(see main text). 


Branch points of h{z) and related qnantities 


n 

Cn 

as 

Rn 

pn 

1 

0.0473642+ i 0.114414 

0.0470313+ i 0.114685 

0.123830 

0.0632582 

2 

0.0177821 + i 0.0673565 

0.0177176+ i 0.067393 

0.0696642 

0.0214698 

3 

0.00946886+ i 0.0475615 

0.0094478+ i 0.0475706 

0.0484949 

0.0114077 


and therefore we dehne 


Cn = c_„ for n < 0 . 


(162) 
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Branch Points 



Figure 7: Branch points Cn’s in the first quadrant of the z-plane obtained as intersections 
of the dashed (red) curve fi with the continuous (blue) curves a 2 n (see appendix\B)- 


Exact values for the first few c„’s are given in the table. An accurate analytic formula (see 
appendix [B]) turns out to be: 


7r(2?T, + 1/2) + iln [e7r(2n + 1/2)] — In [e7r(2n + l/2)]/(7r(2?T, + 1/2))’ ’ 

(163) 

where “In” is the logarithm of real analysis. The relative error with respect to the exact 
value is around 0.3% for n = 1 and quickly goes to zero for increasing n (see Tabled]). 

The hrst derivative of f{w) also vanishes for u —)■ —oo, with w = u + iv and G M. 
That would suggest that also 


z = Q = lim f{u + iv) (164) 

U—^ — OO 

is a branch point of h{z). That however is not true because the origin is not an isolated 
singularity being, as we have just seen, an accumulation point of the order-one branch 
points. 

5.2 The Point at Infinity 

Let us now consider the point at inhnity. We make the usual change of variable ( = 1/z, 
so that 

'll _ 1 _ 77 ; 
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and study the limit C 0. In that limit, u —)■ cxo and 


— (m + iv) 
aw 


(1 - M - iv) - 1 


^ 0 . 


(166) 


The point at inhnity is therefore a branch point of h. Evaluating the higher-order derivatives 
of C{w), one easily hnds that they all vanish for u —)■ oo, implying that = oo is a branch 
point of inhnite order for h{z). Indeed, for u 3> 1, 


f{u iv) 


^u-\-iv _ 2 ^u-\-iv ^u-\-iv{l — l/u) ^u^iv 

_ ^ ^ ~ _ 

u + iv u + iv u u 


(167) 


where the last two relations hold for v u, for example v = 0{1). A sequence of pro¬ 
gressively more crude approximations has been generated above, the last one implying that 
z = oo should be a logarithmic branch point for h. We will see later that h{z) actually 
has an inhnite-order branch point at z = oo, originating from a countable set of order-one 
branch points at inhnity, which ” conspire” to produce a singularity resembling a logarithmic 
branch point. 

Since f{2nin) = 0 for any n ^ 0, /(O) = 1, and f'{2nin) ^ 0 for any n, local inverses 
hniz) do exist because of the inverse function theorem, satisfying 


hn(0) = 27rm for n 7^ 0; ho(l) = 0. 


(168) 


We will analytically continue the above branches to properly cut ^(-planes in the next 
sections. 


5.3 The Branch Hq 

By evaluating both sides of eq. fl258p in 2 ; = 1 (see appendix |^, one hnds that this equation 
relates ho to the principal branch of the Lambert W function: 


hf){z) 



(169) 


ILL has the following series expansion around the origin HI: 


n=l 


n 


n—1 


n\ 


(170) 


which is convergent for 

ici<L ( 171 ) 

The above one is actually the expansion of the bound-state quantum number xid) for 
\g\ -C 1, i.e. for a tightly-bounded particle. Since exp(—l/c„)/c„ = 1/e, by evaluating ho 
in Cn, we obtain 


hoi^Cn) 



(172) 
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Figure 8 : Real plot of the Lambert function w = W{z) for real z > —1/e. Principal 
branch W^: continuous (red) line; zero branch Wq: dashed (blue) line. The two curves are 
separated by the point p = (— 1 /e, — 1 ), where z = — 1 /e is the order-one branch point of 
W and w = —1 its image (see a,ppend,ix{M) . For z 1, W^{z) ~ log(z), while for z < 0, 
1^1 < 1 , Wo{z) K. log(- 2 r). 


The numerical series appearing on the r.h.s. can be exactly summed |14j : 


E 

n=l 


n 


n—1 


n\ e"- 


1 , 


(173) 


implying that 

ho(cn) = 1-^ = dn, (174) 

C-n 

for any n 7 ^ 0. Therefore we conclude that all the order-one branch points of h lie on the 
sheet So, where ho is dehned (see fig J9] for a real plot). Since an order-one branch point 
glues together two different sheets, we have to hnd the second sheet containing the branch 
point Cn for any given n ^ 0. That is accomplished in the next section. 

Another relevant expansion for ho is the one centered in ^ = 1, which is related to the 
case of the ’’loosely bounded particle” discussed before: 

00 

ho{z) = '^kn{z-1)"^, \z-l\<r, (175) 

n=l 

where, by Abel’s theorem, 

r = inf 11 - cJ = 11 - cil = 0.959482. (176) 
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The first few coefficients explicitly read: 


ki 

h 

h 


+ 2 ; 

4 

~3’ 

10 


(177) 

(178) 

(179) 


Function h for real argument 

w = h(z) 



Figure 9: Real plot of the function h for a positive argument, i.e. of the branch w = ho{z) 
for z G M'*'. For z 1, ho{z) ^ log( 2 ;), while for < z ho{z) ^ —l/z. 


5.4 The Branches hn for n 7 ^ 0 

By evaluating a large number of terms (up to 10^) in the power-expansion of hn{z) around 
2 ; = 0 for various n = 1, 2, 3, • • •, 


hn{z) = |z| < R; 

k=0 


and using the standard convergence tests, we find that 


R 


n 


Cn 


1 

_ 

27r|n 


(180) 


(181) 


for any n ^ 0. Let us remark that eq. fll80|) is actually the perturbative expansion of 
in powers of g given in eg. (11261) in a different notation, so we do not repeat the values of 
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the lowest-order coefficients. Since the c„’s are ordered according to their modulus (see 
eq. fllhlD h we conclude that the sheet Si only has the branch point ci (plus an order-one 
branch point at infinity), S 2 has only the branch point C 2 and so on. Therefore: 

hn(Cn) — hiiiCri) = d-a = 1 , Tl 0 (182) 

Cn 


while 


hk{cn) 7^ dn for k ^n. 

The branch point c„ then glues together the sheet Sn to Sq for any n ^ 0. 


(183) 


5.5 Expansions Around the Branch Points c^’s 


In this section we present an expansion of h{z) around the branch point c„ for some fixed 
n ^ 0. Since the latter is an order-one branch point, the series involves semi-integer powers 
of z — Cn- The expansion is obtained by inverting the Taylor series of z = f{w) around dn, 

with Cn = f{dn), 

" - =i; - '-‘’’i"- (184) 

k=2 

We replace on the r.h.s. of the above equation the semi-integer power (Puiseaux) expansion 

OO 

w - dn = 'Y^ {z - Cn)'/^ (185) 

1=1 


expand everything in powers of z — Cn and determine the coefficients in such a way that 
equality holds order by order in z — q 
the explicit expressions: 


where we have defined 


By 




fcS") = 


= 


bi^^ = 


20)^1 . For the first few orders, for example, we have 
1 








2 ’ 


5 (a 






8 ( a. 


An) 


An) 


An) 


du — 


f^Kdn 


k\ 


{k > 2). 


(186) 

(187) 

(188) 

(189) 


we mean an arbitrary but fixed complex square root of for example the 


branch with —vr < arg^; < vr (1^/^ = 1). Changing the convention for \j is equivalent 
to go from one determination of (z — to the other. 


17 


This method is a generalization to the case f'{dn) = 0 of the standard series inversion. 
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To summarize, we have the expansion 

OO 

hnfl{z) =dn + '^ {Z - (190) 

k=l 

The double subscript is to remember that the expansion allows to compute the values 
of both branches hn{z) and for 2 ; sufficiently close to c„. The above series, by a 

generalization of Abel’s theorem, is convergent up to the closest singularity to Cn in the 
z-plane, 

| 2 ;-c„| <p„ = min||c„-c„_i|, |c„-c„+i|| = |c„-c„+i|, for n > 1, (191) 


with Pi = |ci — C 2 I and p_„ 
the table. Asymptotically: 


pn for n < 0. Exact values for the hrst few p„’s are given in 


pn 




1 

27rn{n + 1) ’ 


(192) 


for n S> 1. Let us remark that p„ > 0 for any n ^ 0, but tends to zero for n —)■ ±cxo. 


5.6 Matching Different Expansions 

In this section, in order to check the analytic-geometric structure of the function h previ¬ 
ously established, we match the different expansions obtained. 

Let us denote with D O Cz the convergence region of the expansion of ho{z) coming 
from the power expansion in zero of IT*, 

D = \^z E C\{0} : I exp(—1/2:)/2:| < 1/e}, (193) 

with B (0; i?„) the convergence disk of the expansion of h„(z) around 0, 

Rn) = {z eC:\z\< Rn], (194) 

and by B (c„; p„) the convergence disk of the semi-integer expansion of hnfl{z) in c„ derived 
above (see hgJTOl). 

^ (^ri) pn) — \^Z G C . \z C<fi\ ^ Pn}' (195) 

We can pick up any point p in the double intersection 

pE DnB {Cn, pn) (196) 

and £x the sign of \J~c^ in such a way that, for example, 

^n,o(p) = +ho{p). (197) 

Because of continuity, the sign is constant on any connected component of the intersection 
above. Since the latter is connected, being D and B (cn, Pn) convex, the sign is unique. 
Then for every point 

qE B{0; i?„) n B (c„; p„), (198) 
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Convergence Domains 



Figure 10: Convergence regions of various expansions: 1) disk with continuous (black) 
boundary B{ci,pi): semi-integer expansion for both and hi; 2) “fish-like” region with 
dotted-dashed (red) boundary /3: expansion involving exponentials from the Lambert W 
function for the branch ho; 3) disk with dashed (blue) boundary 5(0, i?i); power expansion 
centered in the origin for the branch hi; f) disk with dotted (green) boundary 5(1, r); power 
expansion of ho^z) centered in z = 1. 


it must hold 

hn,o{q) = -hn{q), (199) 

since similar considerations as those above for the sign also hold in this case. Other checks 
of the analytic structure of h can be obtained by considering the power-expansion in z — 1 
of ho discussed previously, or by taking into account that 

B{Cn]Pn)CB{Cn+l,Pn+l) 7^0, H > 1. ( 200 ) 

Since 5 {cn, Pn) O 5 {cn+ 2 , Pn+ 2 ) = 0, one has to proceed through a chain of disks. In 
general, one can analytically continue the function w = h{z) by moving along paths in the 
z-plane with the ODE derived in the next section or by means of the standard Weierstrass 
procedure involving a sequence of overlapping circles [2T] . 

5.7 Cuts 

We cut the 2 ;-plane of the function h along curves starting from the branch points c„ and 
going to — 00 , of the (arbitrarily chosen) form 

Zn{t) = On - t + i-^^, t>0, (201) 

where c„ = a„ -|- ibn with an^h^ G M, n 7 ^ 0. In fig{TT]we plot the 2 ;-plane cut along fi, 
i.e. the sheet Si where the branch hi is dehned, while in figJT 2 ]we plot the sheet So, where 
the branch ho{z) is defined, having all the cuts n 7 ^ 0. The lifted curves in the 
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Sheet S\ 


0.2 



I—^—,—,—,—I—,—,—,—,—I—,—,—I—^-^—I—I—,—I—,—,—,—,—1 jc 

-0.3 -0.2 -0.1 . 0.1 0.2 


- 0.1 


- 0.2 


Figure 11: Sheet Si of the branch hi, with the branch point ci and the cut ti from ci to 
—oo represented by the dashed (red) curve. 


w-plane are obtained by integrating numerically the ordinary differential equation (ODE) 
for w = h{z), 


dw w 

dz zw — z + V 


( 202 ) 


along tn- A singularity of the solution is expected when the denominator, a polynomial of 
second degree in ^ and w, vanishes |22], i.e. when: 


w = 



(203) 


The above equation, combined with the relation between w and 2 ; provided by eq. (ll52jl . 
gives the equation for the branch points already obtained 


By selecting one determination for the function h, let us define 

Wn{t) = W {Zn{t)) . 

The ODE for T„ explicitly reads: 

dlCji / Wniff) dZn . . 

dt Zn{t) Wn{t) - Zn{t) + I dt 


(204) 


(205) 


(206) 


It is clear that an initial condition to the above ODE for Wn{t) cannot be provided at t = 0: 
solution could not be unique as we would “begin to move” starting from a branch point: 


^n(O) dn, Znift) C, 


(207) 
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Sheet So 


V 



Figure 12: Sheet Sq of the branch Hq, with all the order-one branch points Cn’s of h. To 
construct the Riemann surface ofh, the borders of the cutR, represented by the dashed (red) 
curve going from ci to —oo, have to be glued diagonally with the corresponding borders of 
Si (see the previous figure). All the remaining cuts are plotted as black (continuous) curves 
and have to be connected in a similar way to the associated cuts in the sheets Sn’s. 


As we have just shown, dw/dz, as given by the r.h.s. of eq. (l202p . formally diverges for 
z ^ Cn and w ^ dn- We therefore solve the ODE for t > to with 0 < to 1 0 The initial 
condition Wn(to) can be derived by means of the semi-integer expansion fllQOp 

CXD 

Wn (to) = dn + ^ bj!"^ (Zn(to) - Cn) (208) 

k=l 

For each Zn(to), there are two Wnitofs because of the square root terms, as it should for 
an order-one branch point. The two different initial conditions correspond to the two sides 
of the cut one is moving along to. In practice, to obtain a high accuracy, one has to 
compromise about the numerical value of fg > 0- On fhe one hand, to has to be quite close 
to zero in such a way that Zn{to) — Cn is so small that the series expansion above converges 
quickly, while on the other hand to has to be large enough so that dw/dz is not too big in 
the hrst iterations of the ODE. 


5.8 Asymptotic Expansions 

For I ln„(z)| S> 1, the following asymptotic expansion holds: 

, / X , /X , r, / M Inn [ln„(^)l 1 

hn(z) = \nn(z) -h Ino [ln„(2;)] H-;—+ 


+ 0 


\nn(z) z\nn(z) [lnl{z) 


(209) 


We see here, in a differentiable context, the (algebraic topology) theorem of the unicity of lifted paths, 
once an element in the preimage of a point in the base curve is chosen [23] . 
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Figure 13: Images of the sheets Sn under the branches in the w-plane. hn{Sn), n 0, 
is the region inside the continuous (black) curve Tn, with ho^So) the complementary region 
in the w-plane. Also shown are simple paths allowing one to visit all the Riemann surface 
S ofh. 


where by \n.n{z) we mean the branch of In(^) with argument arg^, in the range 

7r(2?7, — 1) < arg„(2;) < 7r(2?7, + 1), n G Z. (210) 

The dependence on n can be made explicit as: 

Irirfz) = Ino(^) + 27rm. (211) 

For the ’’external” logarithm, we have (conventionally) taken the principal branch lno[- • •]; 
one can actually choose any other branch, as the difference between different choices cancels 
order by order in the expansion, as proved in appendix O An important point is that lnn{z) 
becomes large in two quite different situations: 

1. \z\ 3> F any n. This case is related to the strong-coupling regime of the Winter 
modelO 

2. \n\ S> 1, any z ^ t). This case is related to the high-energy excitations of Winter 
model in the weak-coupling regime, which is our main concern. 

The case h| <C 1 must be discarded because z = 0, as we have seen before, is a non-isolated singularity 
of h, so no expansion around the origin of the above form exists. 
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Actually, in case 1. one can drop the term proportional to l/z, as it is exponentially small 
compared to the other ones, so for \z\ 3> 1, any n, one has: 


hn{z) = \nn{z) + Ino [lnn(2:)] + 


hiQ [ln„(2:)] 

\i\n{z) 


+ 0 


lK(^) 


( 212 ) 


For physics applications (see next section), one needs the asymptotic expansion for a real 
argument. For a; > 0, the expansion explicitly reads: 


hn{,x) = ln(a;) + 27rm + Ino [ln(x) + 27rm] + 

1 


Ino [ln(a;) + 27rm] 
ln(a;) + 27rm 
1 


+ 


X [ln( x) + 27rm] 


O 


In^(^) 


(213) 


where ln(a;) is the real logarithm of a; > 0. As in the general case above, for x S> 1, any n, 
one drops the power-suppressed term proportional to 1/x. The argument 27rn entering the 
above formula can be directly obtained by considering the curve in the ta-plane (see the 
dashed (blue) lines to the right of the imaginary axis in hg JT^ : 


F„ : w = Wniu) = u + 27iin, u > 0. 


(214) 


For n 7 ^ 0, by taking the formal inverse of f{wn{0)) = 0, we obtain h(0) = tan(O) = 27rm, 
implying that we are on the branch n at the initial point of F„: h = hn- As shown previously, 
for M 3> 1, 

f{wM) ZZ (215) 

u 

We now take the inverse again, of the above formula. Since the curve F„ does not cross 
the preimage T„ of the cut in Sn, it is entirely contained in hn{Sn)- By continuity 
therefore the inverse is again hn and u + 2nin ^ hn{e^/u). With x = e^/u, we then 
obtain h„(x) ~ ln(x) -|- 27rm. For n = 0 the argument is similar: one takes the inverse of 
/(wo(0)) = 1 and gets the branch ho; again, the curve Fq does not cross any of the curves 
T„, so it is entirely contained in ho(S'o). For 


X > 0 and n S>-, 

2ttx 


(216) 


the desired expansion reads: 


hn{.x) = 2711 ( n In 


27r 1 n - I X 


O 


In 


— i 


27r(n -I- l/4)x 


2Ti{n -1-1/4) 


+ 


27r(n-hl/4)x ' ^ ( [27r(n-M/4)]' 
For n < 0, one just takes the complex conjugate as: 


(217) 


hn{x) = h_„(x). 


(218) 
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For X < 0, the logarithmic expansion on the r.h.s. of eq. fl209p takes the form: 


hn{x) = In |x| + i7i{2n — 1) + Ino [In |x| + i7i{2n — 1)] + 


O 


Ino [in + i7i{2n — 1)] 
In |x| + m{2n — 1) 


x[ln |x| + m{2n — 1)] 


K(x) 


+ 


(219) 


As in the previons cases, for x -C —1 and n > 1, one drops the power-snppressed term. 
The determination of the argument 7r(2n — 1) in eq. fl219p has been obtained by integrating 
numerically eq. fl202p along the negative axis, i.e. for z = —t with t > 0 with the initial 
condition t(;(0) = 27rm. Let us remark that one cannot evaluate ho(x) for x < 0 by 
integrating the ODE for dw/dz starting from x = 0 because, as we have shown, the origin is 
a non-isolated singularity for this branch and therefore cannot be used as initial condition. 
One cannot even start, for example, from x = 1 (ho(l) = 0), because one cannot pass 
through the origin either. Actually, in the latter case, the following instability occurs. By 
departing arbitrarily small from the real axis, a branch cut tn of some Cn is necessarily hit, 
because the c„’s accumulate at the origin. As a consequence, one leaves the sheet Sq and 
enters the sheet S'„, for some large n, and presumably remains in that sheet, as Sn has c„ 
as its only branch point (apart from the point at inhnity). The consequence is that there 
is a not a well-dehned value of ho(x) for x < 0 as far as ’’physics” (which deals with x G M 
only) is concerned. One has instead to move along great half-circles, starting from positive 
arguments, ending up with a sign depending on the contour orientation. For 


X < 0 and n S> 


1 


27r|x| ’ 


( 220 ) 


the required expansion explicitly reads: 
hn{,x) = 27 ri fn — - j -|-In 


r / 1 \ 1 

In 

1 

to 

1 

'W 

— i —- 


27r(n — l/4)|x| 


27r(n — 1/4) 


+ 


27r(n — 1/4) X 


O 


[27r(n — 1/4)]" 


( 221 ) 


Because of the symmetry between eqs. (12171) and (122ip . one can write both equations in a 
unique form. For 

X e M\{0}, n > (222) 

27r|x| 

it holds: 


hn{x) = 2TTi ( n + 


sign(x) 


-|- In 


27r I n -|- 


sign(x) 


X 


In 


27r(n -|- sign(x)/4)|x| 


27r(n -I- sign(x)/4) 


+ 


27r(n + sign(x)/4) x 


O 


[271 {n + sign(x)/4)] 


2 I ’ 


(223) 


where sign(x) = 1 for x > 0 and —1 otherwise is the function, which returns the sign of x. 
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5.9 Riemann Surface 


By taking into account the results of the previous sections, one finds that the Riemann 
surface S' of h is constructed by gluing diagonally the borders of the cut tn from c„ to —oo 
in Sn to the corresponding ones in S'o, for any n ^ C0. Topologically, S' is a plane, as it is 
usually the case with non-compact Riemann surfaces [2TII2S]. The images of the branches 
^n(S'n) in the tc-plane, the total space of the (branched) covering realized by /, are shown 
in hgJTSl 

One can visit the entire Riemann surface S' by moving along the following simple paths, 
which are also shown in hg{T3l As we have already seen, the dashed (blue) horizontal line 
on the right half-plane of ta, 

r„ : w = Wniu) = u + 27iin, u > 0, (224) 

does not cross any continuous (black) curve T„, passing through dm which is the preimage 
under h of the cut tn- T^ is therefore completely contained in hn{Sn)- Let us remark that 
the point dn divides the curve Tn into two arcs, each of which is mapped by / onto As 
already discussed, these ’’intra-sheet” paths T^ allow one to connect the power expansion 
of hn{z) in z = 0 with the logarithmic expansion at z = cxd. The dashed horizontal (green) 
line on the left half of the w plane, 

Kn : w = Wniu) = u + 27iin, u < 0, n 7 ^ 0, (225) 

crosses the continuous (black) curve T„ once and therefore allows one to go from the sheet 
n to the sheet 0 (or viceversa): 

Kn-. Sn^ So- (226) 

The image curves kn = /(-R„) in the base space z are shown still dashed (in green) in figlTTl 
The dotted-dashed (red) vertical line, 

Jfc : w = Wk{v) = k + iv, u G M, (227) 

where fc is a real constant, crosses all the T^’s for k > Redi, and allows one to go from Sn 
to Sn+i via S'o: 

Jk ■ S'n -H- S'o -H- Sn+ 1 , (n 7 ^ 0; k > Re di). (228) 

For k ^ 1, the image curve in the base space jk = fi-Jk) is basically a great spiral winding 
an inhnite number of times around the point z = 00 , which we have proven to be an 
inhnite-order branch point for h. For fc Re di only some sheets are visited; in this case 
one can still go from any Sn to S'„+i by passing though S'o via the dashed (green) lines. Jo 
coincides with the imaginary axis of the tc-plane, its image in the base space jo = /(^o) 
contains the curves a^’s considered in appendix |B] (see eq. fl27ip i and passes through the 
origin z = 0 an inhnite number of times. Closed curves are obtained in this case and all the 
multivaluedness of h is seen in a very clean way. For n 7 ^ 0, one has indeed fiwoin)) = 0, 
implying /i„(0) = 27rm = rco(’^), while /(tCo(O)) = 1 implies ho(l) = 0 = tCo(O). A portion 
of the image curve jo = /(Jo) in the base space is shown, again dotted-dashed (in red), in 
hgjm 

That is exactly what we do with the function. 

It’s like with the function w = where the origin of the z-plane divides any line passing through 
it into two half-lines, each of which is mapped onto the same half-line of the w-plane. 
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Figure 14: Image curves kn = f{Kn), in the z-plane, of the horizontal dashed (green) lines 
Kn on the left side of the w-plane. By traveling clockwise along kn (which starts in z = 0 
and tends to z = 0 for u —>■ —oo), one goes from Sn to Sq by crossing the cut tn, painted as 
a continuous (black) curve. 


6 Non-Perturbative Analysis 

This is the central section of the paper, in which the limits of perturbation theory are 
rigorously established, non-perturbative results for various observables are obtained and 
compared with the perturbative formulas. The perturbative results for frequencies, widths 
and mixing matrix elements of resonances have all been obtained by inserting into the exact 
expressions the power expansions for k^'^\g) = n — ng , which converge only for 

| 9 | < ( 229 ) 

Therefore, if we consider a fixed resonance (i.e. if we fix n), its properties can be analyzed 
with arbitrary accuracy by means of the perturbative formulas for small enough \g\. Actu¬ 
ally, we are interested in varying n within a specified model, i.e. at fixed g^ by inverting 
the above inequality, the following upper bound on n is obtained: 

Therefore at hxed g ^ 1 (”£xed physics”), only the resonances satisfying the above in¬ 
equality can be described by means of perturbation theory. 

In order to obtain non-perturbative results, we basically insert in the general formulas 
for the observables, the large-n expansion previously obtained, which we repeat here for 


(230) 
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Part of image of the imaginai'y axis of the w-plane, 
Sn 
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Figure 15: Part of the image curve jo = /(Jo); in the z-plane, of the imaginary axis of the 
w-plane, represented by the vertical dotted-dashed (red) line Jq. Traveling counterclockwise 
along the restriction of jo to the interval [ 0 , 27r] (starting in z = 1 and ending in z = 0), 
one goes from So to Si by crossing the cut p. Traveling counterclockwise along /o|[27r,47r] 
(starting and ending at the origin), one goes from Si to So by crossing ti and then from So 
to S 2 by crossing ^ 2 - 


easy of reference: 


hn{-g) = 27ri[n- + in 




4 J 


In 


— z 


27r(n-sign(^)/4)|^| 


27r(n - sign(^)/4) 


+ 


27r(n - sign{g)/4:)g 


O 


[27i{n - sign(^)/4)] 


2 f ’ 


(231) 


where sign( 5 f) is the function which returns the sign of g. The basic equation, which replaces 
the perturbative formula given in eq. fll26p valid for small n, reads: 


^(n)( ) ^ hrrijl 
27ii 


= n 


sign ( 51 ) 


27r 

1 


In 27r|^| 


n 


sign(5() 


In 


\n — sign 


[n — sign(( 7 )/ 4 ] 


471^ g [n - sign( 5 ()/ 4 ] 


O 


[27r(n-sign(^)/4)]" 


+ 


(232) 


The expansion begins with n in both cases bnt, while the pertnrbative expansion is an 
expansion in powers of g, the non-pertnrbative one is basically an expansion in powers of 
1/n. For 

0 < 1 and n 3> 2 T^\g \ ’ (233) 
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the frequencies and the widths take the form: 


(234) 

(235) 



By comparing with the pertnrbative results obtained previonsly, we see that the freqnencies 
have a similar dependence on n to the free case, while the widths grow much less with 
increasing n. 

Since we do not manipnlate perturbative formulas in this section, there is no more 
convenience in introducing renormalization constants and renormalized mixing matrices for 
the resonances, so we switch back to the original (unrenormalized) mixing matrix V{g) and 
resonances On- In terms of the fnnction h, the entries of the mixing matrix V{g) read for 
large n\ 



(236) 


l/2hn(—g) n _j_ 2hn(—g) . 


By inserting the previous expansion for hn{—g), the matrix elements take the explicit form: 






sign(^) 


V{g)i, 


1 + 0 


i n log [27r|5f|?7,] /tt — sign(5f) ?7,/2 — P 


(237) 


A few comments are in order: 1) higher-order corrections in 1/n can be computed in 
a straightforward way; 2) the alternating sign (—1)" is present both in the small-n and 
large-n expansions of V{g)i^n'- h is related to the logarithmic behavior of the fnnction h at 
inhnity. As we have seen in the previous section, this factor is also present in the inverse 
of the pertnrbative expression for V{g), i.e. in Vpt{g)~^. We expect a similar phenomenon 
to occur also in the exact inverse. Indeed, by looking for example at eq. fll48|) . one notices 
that the alternating sign (—1)"' prodnces constructive interference of the harmonics entering 
the sine Fonrier transform aronnd a; = vr, which is the region where non-trivial dynamics 
occurs. For n ^ I with fixed /, i.e. along a row well to the right to the main diagonal. 



n'+> 1. (238) 


Therefore, for large n there is a rather slow power-law decay, with the alternating sign 
(—1)”. The above behavior has to be compared with the pertnrbative one: as we have 
seen, V{g)i^n ~ to 0{g), with convergence properties deteriorating in second order: 
^{ 9 ) 1,11 ~ 1 to (9 (^f^). For I '+> n, i.e. along a colnmn well below the main diagonal: 



( 239 ) 
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A quite slow inverse-power decay with I occurs, with an alternating sign, exactly as in the 
perturbative case. Finally, on the diagonal, one explicitly has: 


^ {.9)n,n 


•i7r/4sign(g) __ 

In {2Ti\g\n) - msign^g)/2 


1 + 0 



I = n. (240) 


Therefore, for large n a square-root divergence basically occurs. In perturbation theory we 
found instead 14,n ~ 1 to 0{g) and 14,n ^ n to O {g"^) (to see that just invert eq. fll25p with 
respect to V {g) or see [6]). 

A numerical analysis (involving the inversion of matrices of various order up to 2 x 10^) 
shows that 


(r(s)-‘) 


l,n 


1 

n 


arg 


(r(s)-‘) 


1,71 + 1 


arg 


(r(9)-‘) 


Ln 


TT for 1. 


(241) 

The hrst relation implies that it is not possible to excite exactly one resonance at a time by 
means of a sensible (i.e. belonging to (R"'")) initial wavefunction and one has necessarily 
to use the approximate (truncated) scheme. Note that 


e{'K-x)Y^ 

n=l 


-sm(na;) 

n 



(242) 


We can summarize the above hndings by saying that an exact treatment of V{g)~^ only 
partially regularizes the perturbative results. The inverse matrix elements {y{g)~^)in 
the non-perturbative case decay faster for n —)■ cxo than in the perturbative case, giving rise 
to less singular initial data. 


7 Conclusions 


The resonances of the Winter model in the weak coupling domain 0 < 4| 1 — which 

is the most interesting one — are subjected to two different regimes. The hrst one is the 
perturbative one, accurately describing the dynamical properties of the n-th resonance up 
to 

The second regime is non-perturbative in character and describes the dynamics of the n-th 
resonance for 

Actually, there is also an intermediate region. 



n ~ 


1 

Sirlsl’ 


(245) 


which, contrary to the above ones, cannot be described by simple analytic formulas. In 
physical terms, we may say that perturbation theory accurately describes the low-energy 
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excitations of the model, while it completely fails to describe the high-energy ones {\g\ -C 1 
always). For the decay widths of resonances, for example, we fonnd: 


~ / 47rc/2n3 + 0{g^) for n < l/(27r|c/|) ; 

\ 2/7rnln (27r|5f|n) +0(1) for n S> l/(27r|5f|). 

The growth of the widths with n is strongly reduced for large n. 


(246) 


—)■ nlnn. 


(247) 


so that the extrapolation of the perturbative formula above in the high-energy region would 
produce completely wrong results. We have also provided explicit non-perturbative expres¬ 
sions for some physically relevant quantities, such as wavevectors, frequencies and mixing- 
matrix entries of resonances, which replace the perturbative formulas previously obtained 
0 E]. We have shown that it is not possible to construct in perturbation theory initial 
wavefunctions exciting exactly one resonance at a time. By inverting numerically the exact 
matrix V{g) truncated up to order N ^ 2x 10^, we have found that that is not even possible 
to construct such a wavefunction in the exact theory. 

The rigorous limitations of perturbation theory, as well as the non-perturbative results 
we found, are all based on the analytic and geometric study of the multivalued function 


w = h{z), 

(248) 

which is the inverse of the entire function 

e"’ - 1 

^ = 

(249) 


where w = 27rik and = —g. It is necessary to complexify the coupling of the model ;s, 
i.e. to leave the physical domain 2 ; G M, in order to go consistently beyond perturbation 
theory. The branch hn{z), n 7 ^ 0, of the function h{z) determines the wavevector of the 
n-th resonance and therefore its frequency, width, etc. A knowledge of all these branches 
is needed to determine the mixing properties of the resonances. The branch ho(^) controls 
instead the properties of the bound state, which occurs in the spectrum for Q < z < 1. 
Each branch hn{z), n 7 ^ 0, is defined on a cut copy of the plane of the complex variable z, 
let’s call it Sn, and has order-one branch points in 


i 


and at infinity. There is a remarkable connection between physics and mathematics: each 
resonance or bound state of the Winter model is associated to a specific sheet of the Riemann 
surface of the multivalued function w = h{z), which is the zero locus of the coefficient b{k, g) 
multiplying the backward wave exp(—ifca;) in the continuous-spectrum eigenfunctions. The 
transcendental equation satisfied by the c„’s reads 



e 


l/z ^ 


e 


(251) 


and has been derived in different ways as: 1 ) the zero locus of the first derivative of the 
inverse function oi w = h{z)] 2 ) the singularity set of the ordinary differential equation for 
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dw/dz] 3) the preimage of the order-one branch point —1/e of the Lambert W fnnction 
nnder the fnnction q{() = — exp{—!/()/(. Qualitative properties of the c^’s have been 
determined, together with explicit (approximate) analytic formulas. The distance of Cn 
from the origin determines the convergence radius of the perturbative expansions (power 
series in z) for the n-th resonance: 

Rn = \Cn\ ~ • (252) 

At variance with respect to all the other sheets, the ’’bound-state” sheet S'o, where the 
branch h^^z) is dehned, contains has all the order-one branch points. The fact that 

—)■ 0 for n —)■ ±oo (253) 

has various physical consequences. The first one is that high-energy resonances (n ^ 1) 
admit a perturbative description for smaller \g\ compared with the lower-energy ones, in 
agreement with physical intuition. The second consequence is that the point = 0 is 
a non-isolated singularity in the sheet S'o, being an accumulation point of branch points. 
Since = 0 is a non-analyticity point, no convergent perturbative expansion for bound-state 
quantities, such as the energy, the wavefunction, etc. exists. We have however provided spe- 
cihc expansions for the bound-state energy: a power expansion m. g + 1, which is convergent 
for \g\ not too small, as well as an expansion for 0 < 1 which involves the functions 

(non analytic in the origin) l/g"^ and e"’/®, n>l{g<Q always). The third and last conse¬ 
quence of eq. (l253p is a severe limitation of perturbation theory: quantities involving in an 
essential way an infinite number of resonances cannot be described by means of a power 
series in g^ the convergence radii RnS of the resonances involved indeed accumulate to zero 
in this case. To excite the /-th resonance, one can use an initial function t = 0; N) 

in which the contributions of the resonances with order n = 1 , 2 , • • • , / — 1, / -|- 1 , • • • have 
been subtracted, with I < N. Let us remark that f = 0; 5 ^, N) can also be computed 

in perturbation theory for \g\ < l/(27riV). 

On the mathematical side, we have constructed the Riemann surface S of the multival¬ 
ued function w = h{z) by gluing together the sheets SnS, where the hn{zys are dehned. 
We found that all the SnS with n 7 ^ 0 are connected via the order-one branch points in Cn 
and at 00 , to Sq. All the resonance sheets Sn, n 7 ^ 0, therefore are not directly connected 
to each other, but only through the bound-state sheet Sq, with the more general coupling 
one could think of, namely square-root branch points. In physical terms, it is remarkable 
that the resonances ’’talk to each other” only indirectly, via the bound-state, which does 
not even appear in the spectrum of the model in the repulsive case, i.e. for g > 0. 
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A The Lambert W Function 

The Lambert W function is dehned as the formal inverse of the entire function 

z = g{w) = w e’”, (254) 

i.e. w = W{z) (see [Hj for a detailed discussion). Symbolically: 

W = g-\ (255) 

The function g is not one-to-one because of the occurrence of the (complex) exponential, 

so W is multivalued. The branch points of W can be determined much in the same way 
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as we have made for the function h. For example, since g'{—l) = 0 while ^ 0, 

the Lambert function has an order-one branch-point in z = g{—l) = —1/e. W also has 
infinite-order branch points in z = 0 and z = cxo, so it can by considered as a generalization 
of the complex logarithm. It is actually one of the simplest multivalued functions having 
both algebraic and logarithmic branch points. There is a branch of W which is analytic in 
the origin, where it vanishes, which is called the principal branch and which we denote by 
VF*. The principal branch is real for a real argument in the range —l/e<z<oo (see 
figJH])- We call Wo the branch connected to IF* via the order-one branch point in z = —Xje. 
The zero branch is real for —1/e < 2 ; < 0. 


A.l Connection to the h Function 


Let us see how IF is related to our function h. Let us assume 2 ; 7 ^ 0 and introduce a variable 
^ = 'C(^) such that: 


By applying h ^ 


h{z) = - 
f on both sides, we obtain: 

ee« = - 


- - e- 

(256) 


g-iF 

(257) 




Since IF = g{^) = ^ e^, the following formal relation holds between the Lambert 
function and the function h\ 


h{z) = 



(258) 


Let us stress however that the above equation involves multivalued functions on both sides 
so, to have an effective utility, one has to specify the branches of both functions. Let us 
also remark that the above equation, as it stands, cannot be directly applied to study a 
neighborhood of 2 ; = 0 of the branches with n 7 ^ 0 , as the latter are analytic in the 
origin. 

Let us see a simple consequence of the above relation. Since the Lambert function 
has an order-one branch point in = — 1 /e, we derive that h{z) has order-one branch points 
whenever 

g-iF ^ 1 
2 : e 

as we have already found directly (see eq. fll56p L 

We can also show that 2 ; = 0 is a non-isolated singularity. Since the function 

q{z) = ^ (260) 

z 

has an essential singularity at the origin, according to Picard’s theorem it basically assumes 
every complex value in any punctured neighborhood of the origin, such as for example 

4(0) = {z e C : 0 < |z| < £}, e>0. (261) 

One can consider paths 7 contained in /^(O) with images ^[ 7 ] going around any number of 
times the branch points of the Lambert IF function. 
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B Evaluation of the Branch Points c^’s of h 


In this section we evaluate the (in general complex) zeroes of the transcendental equation 

exp(—1/z) 1 

z e’ 

which determine the order-one branch points of the function h, with the exception of the 
point z = 1 which, as shown in the main text, is a zero but is not a branch point. 

Let us first present a qualitative argument showing that there cannot be any large c„, 
i.e. any Cn with \cn\ 3> 1, so that all the zeroes fall inside a limited region of the z-plane. 
For \z\ S> 1, one can expand the exponential around the origin to give 



exp(-l/^) 


2 : 



z z"^ 2z^ 


+ 


< 1 , 


(263) 


while the right-hand-side is of order one. Therefore we found a contradiction by assuming 
\cn\ S> 1. The above argument can be refined by taking the imaginary part on both sides 
of ea.fIM: 

exp(—1/^) 


Im 


= 0 . 


By writing z in polar coordinates, 
with 


z = p 


0 < p < cxD and — vr < <p < tt, 

the solution of the above equation is the union of the sequence of curves 


(264) 

(265) 

(266) 


sinp 

• P PniSP) ) 

ip -\r nil 

where n is any integer (see figlTj). Note that 


\pn{p)\ < 


7r(|n| — l) 


—)■ 0 for n —)■ ±oo. 


(267) 


(268) 


implying that the a„’s shrink to zero in that limit. By requiring the radius to be positive, 
we derive the following angular ranges for the a^’s: 


po : (-7r,-f7r); (269) 

Pn : [0, -l-vr) for n > 0 and (— tt, 0] for n < 0. (270) 


Therefore lies in the upper half-plane for n > 0 and in the lower one for n < 0. Since 
the modulus of the r.h.s. of eq. 02671) is < 1, it follows that all the roots of 02621) lie inside 
a closed disk of radius one centered at the origin. Note that 


2 ; = Zn{ip) 


exp(2i(p) — 1 
2i(p + 27rin ’ 


( 271 ) 
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i.e. Un is the image of a segment of the imaginary axis of ta-plane. 

We now prove that the branch-point eqnation (12621) has two infinite seqnences of soln- 
tions, {c,i}n>o and {cn}n<o, both converging to zero and contained in a curve with a vertical 
tangent at the origin. It follows in particular that 


Rec„ 

ImCn 


—)■ 0 for n —)■ ±cx3. 


(272) 


All the branch points lie on the curve, let’s call it /3, obtained by taking the modulus on 
both sides of eq. (12621) . 

(273) 


/5 : cos (p = — p In j 


for 0 < p < 1. As (p = (p(p), let us also set for convenience (p(0) = ±7r/2, where the sign 
is determined by continuity. The curve (3 is symmetric for p —p, has a vertical tangent 
in z = 0, a cusp in z = 1 and bounds a “fish-like” region (see figs 17] and dU]). The last 
condition to impose, to avoid spurious solutions, is 


Re 


exp(—1/z) 


> 0 , 


giving 


, 71 sm p , 71 

27ik -<-p < 27rfc H—, 

2 p ^ 2’ 


for some integer k. By replacing eq. (l267p in the above inequality, we obtain 

2k 

implying that only even n’s are allowed: 


l<n<2k + \, 


(274) 


(275) 


(276) 


n = 2k 


(277) 


with k any integer . By combining the above results, we obtain that any branch point Ck 
lies at the intersection of the curve a 2 k (eq. (l267p l with the curve j3 (eq. (l273p h 

Ck = a2fc n V/c 7^ 0. (278) 

There is an infinite number of such intersections, contained in two sequences converging at 
the origin along /9, as claimed. 


B.l Explicit Solution 

Let us now derive an explicit approximation for c„ for n 3> 1. Since exp(—1/^) has an 
essential singularity in z = 0, it is rapidly varying close to the origin compared to z so, to a 
first approximation, we can fix to an arbitrary value, i.e. set z = C 7^ 0 (the exponential 
never vanishes). We can also assume C = 0(1), as we have already proved that there are 
no large solutions. The dependence on C actually cancels order by order of the iterations 
of the equation derived later in this section. The original equation then simplifies to 

exp(-l/;2:) = -, (279) 

e 
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having the solutions 


(280) 


^n. — 


27rn + ilne/C’ 


with n any nonzero integer. By writing 


C-n — 


27rn + i\\ie/C + kn 

one obtains the (still exact) equation for kn 


_ 


C{27in + i In e/C + fc„) ’ 


( 281 ) 


( 282 ) 


which can be used to set up a recursion scheme generating a series converging to c^. By 
writing 


k — — t(i) 4- t(2) I ^(3) I .., 

1=1 

and assuming (as is verihed a posteriori) 

\n\ > |A:«| > \k^^^\ > \kl^^\ > ••• 
for \n\ S> 1, one obtain for example at hrst order: 


(283) 


(284) 


_ 


C(27rn + i Ine/C') 


(285) 


The second recursion produces the following formula, which is quite accurate and is used 
in the main text (in which we have set at the end C = e)\ 




7r(2?7, + 1/2) + i In [e7r(2n + 1/2)] 


In \eTi{2n + 1/2)] /(7r(2?7, + 1/2)) 


+ 0 


1 

(27™)"^ ’ 

(286) 


for n > 1. For n < 0, one just uses the symmetry relation in eq. 01621) . Note that the 
expansion parameter is actually 1/(2?™), so there is a rather good convergence of the 
above series even for n = 1. The images of the branch points read: 


dn = ivT 



— In 


TT 



In [e7r(2n + 1/2)] ^ 

7r(2n + 1/2) 


1 

(27m)^ 


(287) 


To summarize, we have obtained two sequences of branch points of order one for h for 
n > 0 and n < 0, both converging to zero for n —>■ ±cx) respectively. 


C Absence of Spurious Branches in Asymptotic Ex¬ 
pansions 

Since the function w = h{z) has a logarithmic branch point in z = oo, one expects that the 
expansion around inhnity involves powers of ln(4) only, plus of course eventual single-valued 
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functions of By looking at the asymptotic expansion obtained for h{z), the above fact 
literally is not true, as the composition of the complex logarithm with itself does appear. It 
may seem therefore that the expansion produces spurious multivaluedness. In this appendix 
we show that actually this is not the case, namely that the ’’external logarithm” does not 
give rise to any additional (and nndesired) mnltivaluedness. 

Let us sketch the proof for the simplest non-trivial case, the expansion to first order 
(the hrst two terms). Our expansion. 


hn{z) = ln„(^) + Ino [In,,(2;)] + O 


.ln„(z)J ’ 


(288) 


involving the principal branch of the external logarithm, has to be compared with an 
expansion involving a different branch of the external logarithm; 


hn{z) = \nm{z) + Infc [\nm{z)] + O 


\n„i{z)_ 

where fc is a fixed integer, while \n\ S> 1 by assumption. By setting 

n = m + fc, 

the difference between our formula and the new one reads: 

1 


hn{z) - hn{z) = logo 


1 + 


lno(2;) + 2m{n — k) 


n 


+ 0 -] = 0 - 


n 


(289) 


(290) 


( 291 ) 


i.e. the difference of the logs is absorbed by the higher order terms, as we claimed. 

Let us note that the two expansions actually produce slightly different numerical values. 
In general, there is an ambiguity in snch trnncated expansions, which is never removed, but 
only pushed to higher orders by adding more and more terms. 


That is the generalization to logarithmic branch points of the Puiseaux expansion around a branch 
point a of order n — 1, which involves the integer powers of (z — a)^/" [2D]; actually, it is the formal limit 
for n —>■ 00 . 
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